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ABSTRACT 


_/ 

The  conventional  low-power  CW  doppler  projectile  tracking  radar  is  investigated 
through  simulation.  The  power  requirement  for  full  range  tracking  with  real-time  track 
data  processing  stopped  within  160  m  of  the  point  of  impact  is  obtained.  This  power 
can  be  lowered  bv  a  factor  of  more  than  4  if  a  set  of  three  switchable  filters  is  utilized. 
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I.  INTRODUCTION 


At  a  proving  ground,  it  is  always  desirable  to  record  the  trajectories  of  test-fired 
munitions  for  later  analysis.  This  is  especially  true  when  anomalies  occur.  Continuous 
wave  (CW)  radars  have  often  been  adopted  for  this  purpose  because  of  their  simplicity. 
Compared  to  a  pulsed  radar  of  similar  range  performance,  a  CW  radar  has  no  modulator 
in  its  transmitter,  hence  no  range  gating  is  necessary  in  its  receiver.  Furthermore,  a  CW 
radar  has  a  bandwidth  corresponding  to  the  doppler  shift,  which  rarely  exceeds  100 
kilo-IIertz  (kHz)  while  a  pulsed  radar  has  a  band  width  in  the  mega-Hertz  range.  Thus 
a  CW  radar  can  be  designed  to  carry  out  coherent  signal  integration  over  a  longer  time 
interval.  This  increases  the  integration  gain  and  lowers  the  required  transmitter  power 
for  the  CW  radar. 

High  gain  antennas  of  very  narrow  beamwidths  are  also  used  to  reduce  the  trans¬ 
mitter  power  requirement.  It  is  important  to  aim  such  antennas  accurately  on  the  target. 
To  acquire  the  bearing  and  elevation  angles  of  a  targe’  directly,  several  receiving  feeds 
together  with  complicated  comparator  and  multiple  data  channels  are  needed.  ^Since  the 
trajectory  of  a  test  round  can  be  roughly  predicted,  an  alternative  to  the  actual  acquisi¬ 
tion  of  target  tracking  angles  is  by  steering  the  antenna  according  to  the  predicted  tra¬ 
jectory.  The  directly  measurable  data  will  be  the  doppler  shift  of  the  target  as  a  function 
of  flight  time.  Discrepancies  between  this  data  and  the  doppler  shift  based  on  the  pre¬ 
dicted  trajectory  should  be  carefully  analyzed  and  explained.  When  an  anomaly  occurs 
which  causes  the  target  to  break  off  track,  this  recorded  data  will  be  valuable  for  diag¬ 
nosing  the  problem. 

Many  CW  radar  tracking  systems  are  available.  For  example,  a  model  DR  5S2  X- 
band  Instrumentation  Radar  manufactured  by  TERMA  Elektronik  AS  of  Denmark  is 
being  installed  at  the  Ta-Fu  Proving  Ground  in  Taiwan.  These  systems  typically  track 
a  long-range  projectile  only  up  to  the  summit  of  the  trajectory.  Since  the  point  of  impact 
of  a  long-range  projectile  is  of  great  interest,  it  is  desirable  to  extend  ihe  range  of  the 
radar  to  cover  as  large  a  portion  of  the  trajectory  as  possible  while  keeping  the  system 
relatively  simple.  The  purposes  of  this  thesis  are  to  investigate  the  maximum  range  at¬ 
tainable  using  a  simple  CW  radar  and  to  determine  the  required  minimum  transmitter 


When  tracking  a  projectile,  a  CW  radar  detects  only  the  doppler  frequency  shift  due 
to  the  projectile  motion.  This  frequency  shift  is  proportional  to  the  component  of  the 
projectile  velocity  pointing  radially  away  from  the  radar.  By  combining  this  information 
with  the  tracking  angles,  the  ground  range  of  the  projectile  can  be  deduced  through  in¬ 
tegrating  this  velocity  component.  A  low-pass  filter  can  be  employed  to  perform  this 
signal  integration  in  real  time.  The  output  can  be  used  to  update  the  predicted  track  if 
a  data  filtering  mechanism  such  as  the  Kalman  filter  is  utilized  in  real  time.  Since  in¬ 
formation  on  the  doppler  integrating  filter  is  not  provided  by  any  manufacturer,  one  way 
to  implement  this  low  pass  filter  is  suggested  in  Chapter  2.  The  design  described  therein 
is  consistent  with  the  operational  instruction  for  the  aforementioned  system.  It  is  needed 
for  tracking  simulation  required  in  Chapter  IV.  The  trajectory-updating  Kalman  filter 
is  a  desirable  improvement  to  currently  available  systems  on  the  market.  Its  design  is 
outside  the  scope  of  this  thesis. 

As  explained  above,  trajectory  estimation  is  necessary  for  tracking  with  a  simple  CW 
radar.  Since  the  radar  and  the  gun  are  moving  with  the  spinning  Earth,  the  familiar 
reference  frame  in  which  the  radar  and  the  gun  are  stationary  is  an  accelerating  frame 
which  is  difficult  to  work  with.  Instead,  by  ignoring  the  rotation  of  the  Earth  around 
the  sun.  the  inertial  frame  of  reference  having  its  origin  at  the  center  of  the  Earth  and 
its  z-axis  along  the  spin  axis  of  the  Earth  simplifies  the  equations  of  motion  of  the 
projectile.  The  z-x  plane  can  be  chosen  to  contain  the  radar  at  t  =  0  when  the  projectile 
is  fired.  To  obtain  the  familiar  ground-based  quantities,  transformation  to  the  inertial 
frame  in  which  the  radar  (or  the  gun)  is  passing  instantaneously  through  the  z-x  plane 
can  be  performed.  In  this  frame,  the  distance  and  the  relative  velocity  of  the  projectile 
from  the  radar  (or  the  gun)  can  be  found  easily. 

The  projectile  will  be  spinning  around  its  axis  if  a  rifled  gun  barrel  is  used.  The  lift 
force  due  to  air  pressure  will  rotate  the  projectile  to  align  its  axis  along  its  air  velocity 
vector.  The  kinematics  of  the  projectile,  including  its  angular  momentum,  is  most  easily 
described  in  terms  of  the  motion  of,  and  the  rotation  around,  its  center  of  mass  [Ref. 
1].  The  orientation  of  the  projectile,  which  lies  along  one  of  its  principal  axis  of  inertia, 
can  be  described  in  terms  of  two  of  the  three  Euler  angles  [Ref.  2].  leaving  the  third  angle 
for  measuring  the  spin  of  the  projectile. 

It  can  be  shown  that  the  effect  of  gravitational  force  on  the  projectile  orientation  is 
negligible.  The  ratio  between  the  gravitational  force  on  the  center  of  mass  to  the  force 
affecting  the  angular  momentum  is  of  the  order  of  the  square  of  the  ratio  of  projectile 
size  to  the  Earth  radius  [Ref.  3[.  Thus  the  change  of  the  orientation  of  the  projectile  axis 


is  due  largely  to  the  resistance  of  air.  For  a  streamlined  body,  the  lift  from  the  air  flow 
is  large  when  there  is  an  angle  of  attack  [Ref.  4],  Furthermore,  the  rotation  of  the 
projectile  axis  is  slow.  It  is  reasonable  to  assume  that  the  orientation  of  the  axis  of  the 
projectile  is  always  along  its  air  velocity.  This  assumption  will  be  utilized  to  determine 
the  radar  cross-section  of  the  projectile. 

It  is  easily  seen  in  the  ground  coordinate  that  the  projectile  velocity  is  much  greater 
than  the  wind  velocity  so  that  the  latter  only  introduces  small  corrections  to  the 
projectile  axis.  Since  the  wind  acts  on  the  projectile  over  the  entire  trajectory,  a  better 
estimate  of  the  range  can  be  obtained  by  including  the  wind  velocity  profile  to  trajectory- 
prediction.  A  less  elaborate  alternative  is  to  adopt  a  wind  profile  model  such  as  the 
Ekman  spiral  for  wind  speed  approximation  [Ref.  5],  Since  the  change  in  the  predicted 
range  is  not  expected  to  be  large,  the  wind  velocity  is  treated  as  constant  throughout  this 
thesis. 


3 


II.  PROPOSED  C\\  RADAR  DESIGN 
A.  RADAR  SYSTEM  SCHEMATICS 

The  receiver  design  of  a  radar  is  critical  to  its  range  performance.  Since  no  specific 
information  for  any  C\V  radar  tracking  system  is  available,  the  following  receiver  con¬ 
figuration  is  proposed.  Here  Figure  1  shows  the  schematics  of  the  radar. 


ANTENNA 

BEAM 

INFORMATION 


Figure  1.  Schematics  for  Proposed  CW  Radar  Design. 


The  received  signal  at  the  frequency  from  the  receiving  antenna  is  fed  to  a 

narrow-band,  low-noise  RF  pre-amplifier  with  a  bandwidth  of  1  MHz  centered  at  10.475 
GHz.  The  noise  figure  of  such  a  pre-amplifier  can  be  as  low  as  1.5  to  2  dB,  with  a 
minimum  gain  of  28  dB  [  Ref.  6  j.  This  low-noise  pre-amplifier  can  limit  the  receiver 
noise  figure  (  Fn  )  to  about  2  dB.  The  1  MHz  bandwidth  is  chosen  according  to  the  es¬ 
timated  maximum  doppler  frequency  (fdw)  and  the  fact  that  the  target  is  always  flying 
away  from  the  radar.  The  output  signal  from  the  low-noise  pre-amplifier  is  reduced  to 
the  base  band  at  the  mixer,  and  then  passed  through  the  video  filter  of  bandwidth  Br 
The  filter  output,  at  the  doppler  frequency/,,  can  be  recorded  either  with  an  analog  or 
a  digital  recorder.  The  bandwidth  of  this  recorder  should  be  wider  than  fdmci  to  allow 
truthful  recording  of  the  raw  data  for  later  analysis.  The  doppler  signal  is  then  fed  into 
the  square-law  detector. 

The  output  of  the  square-law  detector  oscillates  at  a  frequency  of  2fd.  Since  this 
frequency  is  proportional  to  the  velocity  of  the  projectile  receding  away  from  the  radar, 
it  can  be  integrated  to  obtain  the  projectile  range  information. 

The  integration  is  most  easily  accomplished  by  counting  the  number  of  cross-over 
points  of  a  pre-determined  level  by  the  square-law  detector  output.  As  to  be  explained 
later,  this  number,  denoted  by  MX,  is  the  displacement  of  the  projectile  away  from  the 
radar  measured  in  half  the  transmitter  wavelengths. 

The  square-law  detector,  a  counter,  and  a  clock  will  be  adequate  for  measuring  the 

range  and  the  velocity  of  the  projectile  receding  from  the  radar.  First,  a  measurement 

base,  MB.  is  set  in  terms  of  the  number  MX  of  transmitter  half  wavelengths.  Then  the 

time  Tn  for  the  output  of  the  square-low  detector  to  run  through  MX  pre-set  level 

V  z> 

crossings  is  measured.  Here  MB  is  the  radial  range  increase  of  the  projectile  and  — - 
is  the  radial  velocity  of  the  projectile. 

As  discussed  in  the  introduction,  the  trajectory  of  the  projectile  is  predicted  before¬ 
hand.  The  elevation  and  bearing  angles  of  the  radar  are  updated  at  the  end  of  every  Tr 
interval.  This  information  is  sent  to  the  servo  control  system  to  keep  the  antenna  on  the 
projectile. 

B.  THEORY 

The  theory  of  operation  for  the  proposed  CW  radar  system  will  be  described.  The 
doppler  effect  is  used  in  a  CW  radar  for  moving  target  detection.  This  effect  was  found 
by  Christian  Doppler  in  1842.  It  states  that  the  frequency  of  a  harmonic  signal  will  shift 
when  reflected  by  a  moving  target.  The  change  in  frequency  is  proportional  to  the  ve- 
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locity  of  the  target.  Let  f,  denote  the  doppler  shift,  Vr  the  radial  velocity  of  the  target 
with  respect  to  the  radar,  c  the  velocity  of  light,  and  f,  the  transmitted  carrier  frequency. 
The  doppler  frequency  shift  is 


a- 2-f 


(2.1) 


2V  — 
r  c 


(2.2) 


The  doppler  is  zero  when  the  trajectory  of  the  projectile  is  perpendicular  to  the  radar 
line  of  sight.  It  is  positive  for  an  approaching  target  (see  Figure  2). 
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Figure  2.  Spectra  of  Received  Signals:  (a)  no  doppler  shift,  no  relative  motion  (b) 
approaching  target  (c)  receding  target. 


The  frequency  of  the  signal  out  of  the  square-law  detector  is  2 fd.  Thus  the  number 
of  crossings,  MX.  of  this  signal  over  a  pre-set  level  within  a  time  interval  T„  is  : 

A/.Y  =  r„  x  (2.3) 

=  7"„  x  2  —r-  (2.4) 

According  to  equation  (2,1),  the  radial  displacement  of  the  projectile  over  this  period 
(Figure  3  on  page  7)  is  : 


i 


6 


Figure  3.  Measurement  Base  and  Velocity  vs.  Time. 


mb  =  r„  x  r ; 

(2.5) 

=  .U.Vxy 

(2.6) 

For  example,  if  the  transmitting  frequency  /  of  interest  is 

10525  MHz,  we  can  set 

the  number  MN  to  100.  Then 

100  x2.907  x10s 

(2.7) 

2  x  10525  x  10* 

=  1.4234  m 

(2.S) 

This  is  a  small  change  in  targe*  location  for  updating  the  radar  pointing  direction. 
Since  MB  is  the  displacement  of  the  target  along  the  radar  pointing  direction,  this  dis- 
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placement  can  be  projected  onto  the  ground.  The  ground  range  can  be  derived  by  add¬ 
ing  up  the  individual  projected  ranges. 
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III.  PROJECTILE  TRAJECTORY  PREDICTION 


To  analyze  the  range  performance  of  the  radar,  information  about  the  trajectory  of 
the  projectile  is  required.  There  are  many  factors  such  as  air  drag,  gravitational  force, 
projectile  shape  and  spin,  Earth  rotation,. ..etc.,  which  can  affect  the  trajectory.  As  dis¬ 
cussed  in  the  INTRODUCTION,  the  projectile  spin  mainly  provides  a  stabilizing 
mechanism  to  the  projectile  and  is  assumed  to  have  little  effect  on  its  trajectory.  The 
direction  of  the  projectile  axis  (or  heading)  is  assumed  to  be  along  the  direction  of  the 
air  velocity  of  the  projectile.  Thus  only  the  motion  of  the  center  of  mass  of  the  projectile 
needs  to  be  considered.  The  projectile  center  of  mass  and  the  center  of  air  pressure  are 
assumed  to  be  the  same.  The  trajectory  is  thus  reduced  to  that  of  a  particle  trajectory. 

Wind  velocity  can  increase,  decrease  or  deviate  the  trajectory  through  air  drag  on 
the  projectile.  In  this  thesis,  it  is  assumed  that  wind  velocity  is  constant  over  the  course 
of  the  projectile;  it  is  usually  called  the  ballistic  wind. 

Drag  force  depends  on  the  air  velocity  of  the  projectile.  It  points  against  the 
projectile  heading.  Its  magnitude  is  proportional  to  the  square  of  projectile  air  speed, 
the  local  air  density,  the  projectile  cross-section,  and  the  drag  coefficient.  Air  density 
depends  on  the  temperature  profile  of  the  atmosphere.  The  drag  coefficient  is  very 
complicated  and,  in  fact,  varies  with  the  projectile  air  speed.  Measured  data  of  the  drag 
coefficient  are  available.  The  drag  force  on  the  in-flight  projectile  is  given  below  : 

7  =  -hj-ACD(V-U')2  (3.1) 

where  IU  is  wind  velocity.  V  is  the  velocity  of  the  projectile.  CD  is  the  drag  coefficient, 
A  is  the  cross-section  area  of  the  projectile,  and  p  is  the  local  air  density.  Here  h  is  the 
heading  of  the  projectile,  which  is  a  unit  vector  in  the  direction  of  V  —  If'. 

The  gravitational  force  is  the  only  other  factor  affecting  the  trajectory  in  an  inertial 
frame  of  reference.  In  the  accelerating  reference  frame  fixed  on  the  surface  of  the  Earth, 
there  is  the  fictitious  Coriolis  force.  It  causes  an  apparent  trajectory  drift.  As  observed 
on  the  ground,  this  drift  can  reach  3.1  km  for  a  trajectory  of  27.4  km  ground  range. 
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A.  THE  TRAJECTORY  IN  A  FIXED  INERTIAL  FRAME 

In  order  to  simplify  the  equation  of  motion,  the  Earth-centered  inertial  frame  (r, 
6,  <p)  is  used.  The  plane  containing  the  radar  and  the  spin  axis  of  the  Earth  when  the 
projectile  is  fired  at  t  =  0  is  the  z-x  plane.  The  radar  is  moving  with  the  Earth  at  an 
angular  velocity  of  Q  .  Its  location  in  this  inertial  frame  is  (r„.  dR,  Clt)  . 

The  components  of  the  projectile  velocity  in  this  reference  frame  are  : 


'=  K 

(3.2) 

re  =  Ve 

(3.3) 

r(  sin  6)<j>  -  I'c 

(3.4) 

The  Lagrangian  formulation  of  the  principle  of  least  action  is  convenient  for  deriv¬ 
ing  the  equations  of  motion  of  the  projectile  [Ref.  7],  The  Lagrangian  of  the  projectile 
in  the  inertial  spherical  coordinate  system  [Ref.  I  ,  8]  (  r,  6,  <f)  ),  without  considering  the 
drag  force,  can  be  expressed  as  : 


L  =  T(r,  6.  d>,  r,  6.  4>)  -  U(r,  6,  0)  (3.5) 

L  =  (»'*  -f  r26 ‘  +  r2<p 2  sin20j  —  mga{  1  — y  )  (3.6) 


where  g  is  the  gravitational  acceleration  at  ground  level, 
equation  of  motion  is  : 

d  cL  _  cL 
di  cr  cr 


For  the  r  component,  the 

(3.7) 


r  =  — g  +  rd2  +r4>2  s'm2d 

For  the  6  component  : 

d  cL  =  cl_ 
di  £6 

— ~(r2d)  =  r24>~  sin  6  cos  6 


(3.8) 


(3.9) 

(3.10) 
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and,  for  the  <f>  component: 


d  cL 
dt  <p 


(3-11) 


sin20)  =  0 


The  nonlinear  drag  force  can  be  included.  The  complete  equations  of  motion  become  : 


d  ,,  ,  a  ,1  , 

-j-K-  -(-)*  + 


,  ,  +  , 

7" )  £  4 - T 


(3.13) 


=  cot  0 1-*  +  r/„ 


(3.14) 


(r  sin  6 1'd)  =  rsinfl/0 


(3.15) 


r  =  a  +  h 


(3.16) 


where  h  is  the  height  from  the  Earth  surface  and  a  is  the  radius  of  the  Earth.  The  ratio 
of  projectile  height  to  the  Earth  radius  can  be  expressed  as : 


J-  =  i  +Jl 

a  1  ‘  a 


(3T7) 


The  projectile  velocity  is  : 


1  —  1  rr  +  1  #6  +  1  cd> 


(3. IS) 


=  rr  +  rdd  +  r  sin  0<j><5 


(3T9) 


and  the  wind  velocity  is  : 


ir  =  ivee+  »’<*»  +  »> 


(3.20) 


Note  that  f,  f.  fe  are  the  components  of  the  drag  force  divided  by  the  mass  of  the 
projectile  : 


(3.21) 
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=  fr  +f0e  +f0i  (3.22) 

The  equations  of  motion  can  be  discretized.  Assume  that  the  gravitational  and  the 
drag  force  acceleration  are  approximately  constant  over  the  period  T„  when  the  projectile 
moves  from  r„  to  r„_„  the  equations  of  motion  are  approximately  : 


rn+\  rn 


=  -  g  —  + 


rnen  +  r„  siW0„</>„ +/„ 


(3.23) 


2  A  2A 

T„ 


=  r\  sin  6n  cos  6n<t>2„  +  rj^n 


(3.24) 


rn+\  sin:0„+1d>„+1  -  r2n  sm28n(bn  .  ^  „ 

- T -  =  rn  Sm  0Jtn 

1  n 

From  the  above  equations,  if  rn,  6„.  <t> ,  and  r„.  0„.  <f>r  are  known,  then  r„_,, 
r„.i,  0„.l,  and  0„_,  can  be  derived  as  follows  : 

With  measured  T„,  equation  (1)  provides  ,  then 


(3-25) 

(Li*  </>„-!• 


rn-t] 


—  rn  + 


>'n  +  'Wl 


Equations  (3.26)  and  (3.24)  give  0„.,,  then 


(3.26) 


a  -  e  +  +  j 

Wn+1  ~  +  2  n 


Equations  (3.27)  and  (3.23)  give  then 


(3.27) 


<t>n+ 1  ■  6,+  *"*/"**  fn  (3.28) 

Note  that  the  muzzle  location  is  re  and  the  muzzle  velocity  of  the  projectile,  which 
is  assumed  to  be  known,  is  I',.  Through  the  above  procedures,  the  full  trajectory  of  the 
projectile  can  be  constructed  by  computer  simulation. 
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B.  GROUND-OBSERVED  TRAJECTORY 

In  order  to  obtain  the  trajectory  as  observed  from  the  ground,  the  position  of  the 
projectile  should  be  converted  into  its  relative  location  from  the  radar  antenna  and  from 
the  gun.  The  necessary  coordinate  transformation  matrices  are  given  below  (Figure  4 
on  page  13). 


Figure  4.  Inertial  Coordinate  Systems. 

1.  Trajectory  Relative  to  Radar 

To  set  the  elevation  and  bearing  angles  for  the  radar  and  to  estimate  the  doppler 
shift,  the  location  and  velocity  of  the  projectile  relative  to  the  radar  as  viewed  from  the 
ground  are  needed.  On  the  ground,  the  projectile  motion  is  most  easily  described  in  an 
inertial  frame  in  which,  at  time,  t,  the  radar  is  passing  through  the  z'  —  x'  plane  instan¬ 
taneously.  The  z'  axis  is  chosen  so  that  it  coincides  with  the  z  axis.  The  transformation 
from  the  fixed  inertial  spherical  coordinate  system  to  this  (jr',  y'  ,z')  system  can  be  ac¬ 
complished  through  the  following  rotation  matrices  R ,,  Rrtc ,  and  Rc>, : 
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R,m  = 


(3-29) 


cos  ill  sin  Hr  0 
-  sin  Qt  cos  Q?  0 
0  0  1 


Kec(<t>p) 


cos  <j)p  0  —  sin  4>p 
sin  4>p  0  cos  4>p 

0  1  0 


RcytOp)  « 


sin  6p 
cos  ep 
0 


cos  dp  0 
—  sin  dp  0 

0  1 


(3.30) 


(3.31) 


a.  Velocity  Transformation 

The  projectile  velocity  in  the  (x',y\  z')  coordinate  system  can  be  obtained 
by  first  changing  the  results  in  section  II 1A  from  the  spherical  coordinate  system  into 
the  rectangular  coordinate  system  (x.  y.  z)  : 


rp* 

V 

»»■ 

=  Rrec(4>p)R(>fdp) 

_  f  P°. 

(3.32) 


Since  the  (,v',y,  :')  coordinate  system  is  obtained  by  rotating  the  (x.  y,  z) 
coordinate  system  by  an  angle  ill  around  the  z-axis. 


VP* 

VPX 

ypy 

=  R:m 

lAv 

Vpr 

VPZ 

-  Rrec&p-WRM 


Pr 


V 


V, 


pc 


The  velocity  of  the  radar,  I  in  the  (x',y'  coordinate  is  : 


(3.34) 
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0 

fV 

~  RretWWtyft id 

0 

_  *  Rf  _ 

sin  6  R 

cos  6^ 

o' 

0 

= 

0 

0 

1 

0 

cos  Br 

—  sin  0R 

0 

*  R<f> 

0 

=  1  Ro 

0 

Hence 

Vr  -  V**' 

=  Q(rRsineR)y' 

b.  Projectile  Location  Relative  to  the  Radar 

The  relationship  between  (.t'.ji',  z')  and  (r,  0,  4>)  is  given  by  : 


(3.35) 


(3.36) 


(3.37) 


(3.38) 


(3.39) 


X 

r  sin  6  cos  <5 

y 

= 

r  sin  0  sin  <j5 

_  z  _ 

r  cos  6 

V 

}■' 

=  RJQr) 

X 

V 

2 

-  -  - 

r  sin  6  cos (4>  -  ft/) 
r  sin  6  sin( <i>  -  ft/) 
r  cos  6 


(3.41) 


(3.42) 
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c.  Radar  Coordinate 


Substitute  ( rR ,  dR.  Qr)  into  equation  (3.42). 

rR  sin  0R 
0 

1 R  COS 


(3.43) 


2.  Trajectory  Relative  to  Gun 

Relative  to  the  gun  (muzzle),  a  similar  inertial  frame  (.v  "  .  v  "  , z  "  )  in  which  the 
gun  is  instantaneously  passing  through  the  z  "  -  x  "  plane  at  time  t  can  be  used.  The 
c  "  axis  coincides  with  the  z-axis.  The  coordinate  transformation  is  : 

r  sin  0  cos(0  —  4>c  —  Qi) 
r  sin  6  sin(<5  —  <?c  —  Qt) 
r  cos  6 


The  firing  plane  of  the  gun  is  the  plane  containing  its  barrel  and  passing  through 
the  center  of  the  Earth.  Tor  artillery,  it  is  customary  to  use  the  coordinate  system 
(x,.v,. z;)  in  which  the  gun  muzzle  is  the  origin  and  the  z,  —  xf  plane  is  the  firing  plane. 
Here  z,  is  the  vertical  direction  and  .t.  is  along  the  horizontal  from  the  gun.  The  relation 
between  (.v..  r.,  z.)  and  (,v  "  .y  "  .  r  ")  depends  on  the  gun  bearing  angle  which  is  the 
angle  of  the  .t.  axis  measured  counterclockwise  from  the  East.  If  the  the  gun  is  oriented 
along  6aZC  and  is  located  at  rc  ,  60.  4>c,  the  vector  rF  —  rG  in  the  coordinate  (w.y..  zA  is  : 


0  cos  0 A:C  sin  6a:Q  cos  0  sin  0C  x  " 
0  —  sin  6A:G  cos  6 a,g  x  0  1  0  y  " 

10  0  —  sin  8g  0  cos  0G  z  " 


(3.45) 


Note  that  rt  -  f0  is  the  relative  location  of  the  projectile  from  the  gun.  Hence 
its  components  can  replace  (x  "  ,y  "  ,z  "  )  directly  in  equation  (3.45). 
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IV.  RANGE  ANALYSIS  AND  POWER  REQUIREMENT 


A.  RANGE  EQUATION 

The  power  requirement  of  the  CW  radar  can  be  obtained  from  the  radar  equation 
[Ref.  9]  : 


P ,  = 


(Ai:)hTsrnL(-Y 

G2o/.:r 


(4.1) 


In  equation  (4.1),  the  radar  cross-section  (a),  the  maximum  range  (R),  and  the 
propagation  factor  (F)  will  be  discussed  in  this  chapter  in  detail.  The  Boltzman  constant 
k  is  287.05  n-m  kg-°A.  The  absolute  temperature  T,  is  290CA’  .  The  wavelength  /.  is 
0.02S5m.  Other  parameters  are  based  on  the  configuration  of  the  CW  radar  as  proposed 
in  Chapter  2.  They  are  defined  as  follows. 

1.  The  Receiver  Noise  Figure 

The  receiver  noise  figure  (A)  is  the  ratio  of  the  signal-to-noise  power  ratio  at  the 
receiver  input  to  the  signal-to-noisc  power  ratio  at  the  receiver  output,  assuming  that  the 
receiver  temperature  is  290‘A  .  The  F„  value  for  the  receiver  considered  is  assumed  to 
be  2  dB. 

2.  System  Loss  Factor 

The  loss  factor  consists  of  the  fractional  loss  of  transmitter  power  delivered  to 
the  transmitting  antenna  (A,),  the  fractional  loss  of  the  receiving  antenna  output  power 
delivered  to  the  receiver  input  (A,),  and  miscellaneous  loss  (A„).  They  are  discussed  be¬ 
low. 

a.  Transmission  Loss 

The  transmission  loss  is  the  loss  in  the  transmitting  or  the  receiving  portion 
of  the  transmission  line.  Since  the  transmitter  is  light-weight  and  is  to  be  mounted  very 
close  to  the  antenna,  A,  is  estimated  to  be  about  1.2  dB  due  to  the  use  of  narrow  band¬ 
width  rotary  joints.  The  receiver  pre-amplifier  can  be  mounted  directly  to  the  receiving 
antenna  outpuy  followed  with  the  mixer  immediately.  This  arrangement  will  reduce  A, 
to  about  the  same  level  as  A.. 
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b.  Miscellaneous  Loss 

The  miscellaneous  loss  (Lm)  for  this  CW  radar  system  is  mainly  due  to  the 
quantization  loss  at  the  comparator.  Since  only  two-level  (crossing)  comparison  is 
made,  it  is  not  expected  to  be  significant  and  will  be  ignored. 

Base  on  the  above  argument,  the  total  loss  is  estimated  at  2.4  dB. 

3.  Signal-to-Noise  Ratio 

For  a  tracking  radar,  it  is  important  that  the  receiver  output  signal  is  sufficiently 
high  above  the  noise  level  so  that  information  carried  by  the  signal  can  be  extracted  re¬ 
liably.  With  the  proposed  receiver  configuration,  the  signal-to-noise  ratio  (S  N)  of  the 
output  from  the  video  filter  is  the  figure  of  significance.  First  of  all,  this  signal  is  re¬ 
corded  for  later  analysis.  Its  minimum  S  N  has  to  be  assured.  In  real-time,  after  passing 
through  the  square-law  detector,  this  signal  generates  the  timing  sequence  for  updating 
the  radar  pointing  direction  and  for  evaluating  the  projectile  radial  velocity.  It  drives  the 
gain  control  to  keep  the  reference  signal  at  a  preset  level.  It  is  also  the  input  into  the 
comparator  following  which  the  number  of  signal  cross-overs  is  tallied  and  the  elapsed 
time  is  clocked  as  explained  in  section  IIA. 

The  square-law  detector  and  the  comparator  perform  threshold  detection.  The 
reference  signal  acts  as  the  threshold.  The  comparator  outputs  a  +  1  signal  when  the 
detector  output  is  above  the  reference  level;  it  outputs  a  -1  signal  when  the  detector 
output  falls  below  the  reference.  The  reference  signal  level,  which  is  set  at  a  fraction  of 
the  detector  output  signal,  determines  the  accuracy  of  the  real-time  output  of  the  system. 
In  what  follows,  the  minimum  S  N  ratio  and  the  threshold  setting  in  terms  of  the 
threshold-to-signal  power  ratio  are  discussed. 

The  signal  input  to  the  comparator  is  oscillating  at  twice  the  dopplcr  frequency. 
Within  each  cycle  of  the  signal,  the  portion  rising  above  the  threshold  level  to  the  peak 
and  back  to  the  threshold  level  can  be  considered  as  a  pulse.  This  pulse  is  contaminated 
with  noise,  which  is  assumed  to  be  white.  The  probability  of  each  contaminated  pulse 
falling  under  a  certain  threshold  level  is  well-known.  This  probability  of  missed  de¬ 
tection  is  the  fraction  of  signal  pulses  failed  to  be  detected  by  the  comparator.  Between 
two  pulses,  the  number  of  false  alarms  is  the  number  of  extra  cycles  counted.  These  two 
effects  tend  to  cancel  each  other  so  the  optimum  threshold  setting  is  achieved  by  equal¬ 
izing  the  fraction  of  missed  pulse  and  the  fraction  of  noise  pulses.  The  program  to 
compute  these  fractions  are  listed  in  Appendix  D  [Ref.  10). 

Figure  5  on  page  19,  Figure  6  on  page  20,  Figure  7  on  page  20.  Figure  8  on 
page  21,  and  Figure  9  on  page  21  show  the  fractional  errors  due  to  missed  pulses  and 
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noise  pulses  respectively.  If  the  signal-to-noise  ratio  is  10  dB,  the  comparator  output 
error  is  5.55%  when  the  threshold-to-signal  power  ratio  is  set  at  0.455.  Stronger  signal 
allows  the  threshold  to  be  set  at  a  lower  fractional  value  with  an  improved  error  ratio. 
For  example,  for  a  signal  12  dB  above  the  noise,  the  threshold  setting  of  0.38,  or  -4.2 
dB  below  the  signal,  will  limit  the  error  to  1.15%.  In  this  thesis,  the  minimum  signal- 
to-noise  ratio  is  chosen  to  be  10  dB. 
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Missed  Detection  (+)  &  Noise  Error  (o) 


Figure  6.  Missed  Detection  and  Noise  Error  vs.  Threshold  to  Signal  Power  Ratio 
S  N=  10.5  dB. 


Figure  7.  Missed  Detection  and  Noise  Error  vs.  Threshold  to  Signal  Power  Ratio 
S  \  =  11  dB. 
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Missed  Detection  (  +  )  &  Noise  Error  (o)  *  Missed  Detection  (+  )  &  Noise  Error  fo) 


0.36  0.38  0.4  0.42  0.44 

Threshold-to-Signai  Powe 


8.  Missed  Detection  and  Noise  Error  vs.  ThreshoId-to-Signal  Power  Ratio 
S  N=  11.5  dB. 


Figure  9.  Missed  Detection  and  Noise  Error  vs.  Threshold-to-Signal  Power  Ratio 
S  N=  12  dB. 


4.  Transmitting  and  Receiving  Antenna  Gain 

Separate  but  identical  transmitting  and  receiving  antennas  are  used  to  reduce 
direct  coupling  of  the  transmitter  power  into  the  receiver.  Since  this  separation  is  small, 
it  is  neglected  in  the  analysis.  In  x-band,  antennas  with  one-degree  beamwidth  are  easily 
achievable.  Based  on  this  beamwidth.  the  antenna  gain  (G)  [Ref.  9]  is  : 


^  27000 

G= 

(4.2) 

=  AA.ldB 

(4.3) 

and  the  diameter  D  of  such  an  antenna  is 

c-72*ir 

(4.4) 

-  72  x  *-nf 

(4.5) 

=  2m 

(4.6) 

5.  Video  Bandwidth  and  Filters 

The  required  minimum  video  bandwidth  (Z?,)  is  a  function  of  the  doppler  fre¬ 
quency.  As  previous  mentioned,  the  bandwidth  should  be  greater  than  fame,  .  Although 
the  doppler  frequency  fa  of  the  projectile  increases  rapidly  right  after  firing,  it  decreases 
monotonically  afterwards  by  a  factor  of  more  than  4  over  its  trajectory  .  An  example  is 
shown  in  Figure  10  on  page  23. 
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x104 


Figure  10.  Doppler  Frequency  vs.  Slant  Range. 

Referring  to  the  range  equation  (4.1),  in  order  to  use  the  lowest  transmitted 
power  possible  for  full  range  tracking,  the  radar  range  and  video  bandwidth  can  be 
traded  off.  Note  that,  assuming  a  constant  radar  cross-section,  the  transmitter  power 
depends  on  the  product  A  filter  of  narrower  bandwidth  may  be  employed  when 

the  projectile  reaches  a  greater  distance.  Adoption  of  a  set  of  three  switchable  doppler 
filters  is  recommended.  The  following  is  an  example  for  its  implementation  based  on  the 
doppler  frequency  of  the  projectile  given  in  Figure  10. 

The  maximum  doppler  shift  fdmo,  from  the  projectile  shown  in  Figure  10  is  a 
little  less  than  54  kHz.  The  maximum  range  R  is  about  27  km.  Three  filters  are  to  be 
used  over  three  different  ranges.  Each  of  the  filters  will  have  a  bandwidth  of  1.2  times 
the  maximum  expected  doppler  frequency  over  the  range  in  which  it  is  to  be  used  to  al¬ 
low  for  a  20%  estimation  error.  The  transmitter  power  level  will  be  minimized  if  the 
fourth  power  of  the  maximum  range  times  the  maximum  doppler  shift  in  the  corre¬ 
sponding  range  is  a  constant. 
a.  First  Filter 

When  the  projectile  is  within  16.06  km  slant  range,  the  video  bandwidth  is 
chosen  to  be  1.2  x  fdnci.  Thus  : 
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Bv  =  1.2  xfdmax 

(4.7) 

=  63.9  kHz 

(4.8) 

b.  Second  Filter 

When  the  projectile  moves  beyond  16.06  km.  but  within  23.73  km  slant 
range,  the  doppler  frequency  has  a  maximum  value  of  about  22.03  kHz.  The  video 
bandwidth  is  : 

=  1.2  x  22.03  kHz  (4.9) 

=  26.44  kHz  (4.10) 

c.  Third  Filter 

When  the  projectile  reaches  23.73  km  slant  range,  the  maximum  doppler 
frequency  is  about  12.7  kHz.  The  video  bandwidth  is  : 

Bv  =  1.2  x  12.7  kHz  (4.11) 

=  15.2  kHz  (4.12) 

Since  R%„e,  is  now  less  than  1  4  of  the  original  value  when  only  one  video  filter  is  used, 
the  minimum  required  transmitter  power  will  be  reduced  by  more  than  a  factor  of  4. 
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B.  TRAJECTORY  SIMULATION  AND  RESULTS 

Based  on  the  procedures  for  trajectory  estimation  presented  in  Chapter  III.  a  simu¬ 
lation  program  was  written  and  is  included  as  Appendix  A.  It  was  written  in  PC-Matlab 
3S6.  The  program  is  described  as  follows. 

1.  Input  and  Parameters. 

a.  Earth 

The  operation  is  assumed  to  take  place  in  Taiwan.  The  radar  is  located  at 
the  longitude  of  1 2 1 c 36'27  "  East  and  the  latitude  of  24°36'3S  "  North.  The  Earth  ra¬ 
dius  is  approximately  6.3677225  x  10' m  (average  value  of  the  artic  radius 
6.356912  x  106  m  and  the  equator  radius  6.37S533  x  106  m).  The  gravitational  acceler¬ 
ation  at  sea  level  is  9. SI  m/ sec-’.  The  angular  velocity  of  Earth  spin  is 
7.2921 1S44  x  10  -' rad  sec  [Ref.  11]. 

b.  Radar 

The  transmitting  frequency  is  10525  MHz.  The  measurement  base  MB  can 
be  varied  to  change  the  update  rate  of  the  antenna  pointing  direction.  For  both  the  re¬ 
ceiving  and  the  transmitting  antenna,  the  beamwidth  is  lc  and  the  gain  is  44.3  dB.  The 
antenna  height  and  the  distance  between  the  gun  and  the  radar  can  affect  the  initial 
readings  of  the  do;  pier  signal. 

c.  Gun 

The  firing  elevation  angle  is  assumed  to  be  45c  and  the  azimuth  is  counter¬ 
clockwise  from  the  East  and  is  also  assumed  to  be  45°.  Both  the  elevation  angle  and  the 
azimuth  angle  can  be  varied.  The  distance  between  the  gun  and  the  radar  can  be  ad¬ 
justed.  (Variations  of  up  to  2000  m  are  possible  at  the  test  site.)  The  barrel  length  of 
the  gun  under  consideration  is  7  m. 

d.  Projectile 

The  target  is  a  155  mm  calibre  artillery  projectile.  Its  weight  is  45.4  kg,  with 
a  total  length  of  938  mm.  The  maximum  muzzle  velocity  is  850  m  sec. 

c.  Environment 

( 1  >  Temperature  and  air  density  profiles.  The  temperature  T„  at  sea  sur¬ 
face  is  assumed  to  be  28S°A  and  the  temperature  at  the  projectile  height  h  is  based  on 
the  following  profile  [Ref.  12]  : 


Th  =  Tc  -  a  f 


(4.13) 


where  h  is  the  height  from  sea  level,  aT  is  the  temperature  lapse  rate  (-0.0065  °K  m). 
The  temperature  (Th)  is  used  for  calculating  the  air  density  (ph)  at  the  projectile: 


Ph  =  Poi-T-)1’1^**'2 


(4.14) 


In  equation  (4.14)  p„  is  the  air  density  at  sea  level  (1.225  kg  m3),  and  R  is  the  gas  con¬ 
stant  (2S7.05  n-m  kg-°A'). 

f2j  Wind  profile.  The  wind  is  assumed  to  be  constant  during  the  flight 
time  over  the  projectile  trajectory.  The  value  used  is  the  so-called  ballistic  wind.  It  can 
be  expressed  as  : 

w  =  wr?  +  wee  +  wj  (4.15) 


In  this  thesis,  the  wind  is  6  m  sec,  North-West.  It  should  be  noted  that  the  extent  to 
which  the  estimated  trajectory  can  be  modified  by  the  actual  wind  profile  is  subject  to 
further  study. 

Mach  number.  The  Mach  number  is  the  ratio  of  the  projectile  air 
speed  to  local  sound  speed.  The  sound  speed  at  the  height  of  the  projectile  is  : 

ah  =  N  KR 1).  (4.16) 


where  K  is  1.4  for  idea  gas.  The  Mach  number  (M)  is  given  by  : 


M  = 


i~-  in 


(4.  IS) 


where  I  V  —  IF  |  is  the  projectile  speed  with  respect  to  air. 

/  Drag  Coefficient 

The  drag  coefficient  (CD)  is  required  for  calculating  the  drag  force.  It  de¬ 
pends  on  the  shape  of  the  projectile.  The  Q,  curve  is  measured  experimentally.  It  is 
given  as  a  function  of  Mach  number  M.  For  a  155mm  projectile,  the  CD  curve  is  shown 
in  Figure  1 1  on  page  27. 

2.  Results  of  Simulation 

The  results  of  trajectory  simulation  are  described  below.  The  variable  names 
used  in  the  program  are  given  in  the  parentheses  when  they  are  first  mentioned. 
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a.  Radar  Observed  Trajectory 

The  slant  range  (SRR)  is  the  distance  from  the  radar  antenna  to  the  in-flight 
projectile.  The  curve  showing  the  projectile  height  (Height)  versus  its  slant  range  from 
the  radar  is  shown  in  Figure  12  on  page  2S.  The  total  flight  time  is  87.4  seconds,  and 
the  maximum  range  is  27.4  km. 

b.  Gun  Based  Trajectory 

The  ground  range  is  measured  from  the  gun  to  the  projectile.  The  projectile 
height  vs.  in-flight  time  (t)  and  the  projectile  height  vs.  ground  range  are  shown  in 
Figure  13  on  page  28  and  Figure  14  on  page  29.  The  maximum  ground  range  (Grange) 
is  27.4  km. 


Figure  11.  Drag  Coefficient  vs.  Mach  Number  for  155  mm  Projectile. 


projectile  height 


Figure  12.  Projectile  Trajectory:  Height  vs.  Slant  Range. 


Figure  13.  Projectile  Trajectory  :  Height  vs.  Time. 
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Figure  14.  Projectile  Trajectory:  Height  vs.  Ground  Range. 


The  cross  range  (CRangc)  measures  the  drift  of  the  projectile  from  the  firing 
plane.  The  firing  plane  is  the  vertical  plane  containing  the  gun  barrel.  In  the  case  con¬ 
sidered,  the  maximum  cross  range  is  about  6.3  km  South-East  of  the  firing  plane.  The 
point  of  impact  is  about  3.14  km  away  from  the  firing  plane.  The  curve  of  cross  range 
vs.  ground  range  is  shown  in  Figure  15  on  page  30. 


Figure  15.  Projectile  Trajectory:  Cross  Range  vs.  Ground  Range. 
c.  Radar  Tracking  Parameters 

The  azimuth  (traz)  and  elevation  (trel)  antenna  pointing  angles  are  required 
for  tracking.  The  curves  of  these  angles  vs.  in-flight  time  when  the  gun-radar  distance 
is  500  m  are  shown  in  Figure  16  on  page  31  and  Figure  17  on  page  31. 

The  azimuth  and  elevation  angular  speeds  of  the  antenna  are  shown  in 
Figure  18  on  page  32.  and  Figure  19  on  page  32. 


Figure  18.  Radar  Azimuthal  Tracking  Angular  Speed  vs.  Time. 


Figure  19.  Radar  Elevation  Tracking  Angular  Speed  vs.  Time. 
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The  antenna  servo  system  has  to  point  the  antenna  to  within  half  a  beamwidth  in  the 
direction  of  the  projectile.  Thus  the  half  beamwidth  divided  by  the  angular  speed  of  the 
antenna  is  the  required  response  time  of  the  servo  system.  The  servo  response  time  vs. 
the  projectile  in-flight  time  is  shown  in  Figure  20  on  page  33.  From  the  above  results, 
it  can  be  seen  that  over  the  first  few  seconds,  the  elevation  angular  speed  is  very  high. 


Figure  20.  Servo  System  Response  Time  vs.  Time. 


d.  Doppler  Frequency  and  Velocity 

The  doppler  frequency  {fj  is  a  directly  measurable  data  by  a  C\Y  radar. 
The  doppler  cycle  can  be  detected  with  the  square  law  detector.  In  the  proposed  CW 
radar  configuration  of  this  thesis,  the  doppler  cycle  is  derived  from  the  time  interval  Tn 
for  each  measurement  base  traveled.  This  timing  sequence  has  to  be  simulated  in  the 
tracking  program.  In  the  simulation  program,  the  doppler  velocity  (1^,.)  is  calculated 
from  the  projectile  velocity  relative  to  the  radar  and  the  radar  pointing  direction.  Then 
T„  is  estimated  from  the  time  required  to  travel  the  distance  MB  with  this  doppler  ve¬ 
locity.  The  curves  of  doppler  frequency  vs.  slant  range,  doppler  frequency  vs.  time,  and 
doppler  velocity  vs.  time  are  shown  in  Figure  10  on  page  23,  Figure  22  on  page  34,  and 
Figure  21  on  page  34. 
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Doppler  Frequency  (Hz) 


Based  on  the  above  results  from  trajectory  simulation,  the  projectile  radar 
cross  section  and  the  transmitter  power  can  be  estimated.  They  will  be  discussed  later. 


x104 


Figure  21.  Doppler  Frequency  vs.  Time. 


Figure  22.  Doppler  Velocity  vs.  Time. 


C.  PROJECTILE  RADAR  CROSS-SECTION 

The  projectile  radar  cross-section  (RCS)  is  denoted  by  a.  It  is  the  area  intercepting 
the  amount  of  plane  wave  power  which,  when  scattered  equally  in  all  directions, 
produces  an  echo  at  the  radar.  In  other  words,  the  RCS  is  given  by  : 

a  =  -4 rr /• 2 (  ~r  )2  (4. IS) 

where  £,  is  the  scattered  electric  field  strength  from  the  target,  r  is  the  distance  from  the 
target  in  its  far  field  where  the  quantity  £  is  measured,  and  £,  is  the  incident  plane  wave 
field  strength  at  the  target. 

Determination  of  the  RCS  of  a  complex  target  is  difficult.  Usually  the  radar  cross- 
section  of  a  target  depends  on  the  following  factors:  a)  target  size,  b)  target  shape,  c) 
radar  frequency,  d)  target  aspect  angle,  and  e)  radar  polarization.  The  RCS  is  greatly 
simplified  if  the  target  size  is  large  compared  to  the  radar  wavelength.  Consider  the  ra¬ 
dar  cross-section  of  a  simple  sphere.  The  RCS  as  a  function  of  the  circumference  of  the 
sphere  measured  in  wavelengths  (2 nal/.)  is  shown  in  Figure  23  on  page  36.  The  region 
where  the  target  size  is  large  compared  to  the  radar  wavelength,  i.e.,  where  2 nal/.  is 
greater  than  10.  is  called  the  optical  region  (Ref.  9].  The  RCS  in  this  region  is  roughly 
the  physical  cross-section  at  area  of  the  target. 

The  155  mm  projectile  is  in  the  optical  region  at  the  transmitter  frequency 
(10525  MHz)  of  interest.  Its  RCS  value  will  be  approximated  with  its  geometrical 
cross-sectional  area.  Because  of  the  rotational  symmetry  of  the  projectile,  the  radar 
cross-section  of  the  projectile  is  a  function  of  the  angle  Q,..  between  the  radar  pointing 
direction  ( D )  and  the  projectile  heading  (/? ). 

The  radar  pointing  direction  D  is  given  by  : 

D  ~  -iLZJj—  (4T9) 

I  >P ~  >R  I 

where  rr  is  the  projectile  position  vector,  rR  is  the  radar  position  vector  as  defined  in 
Chapter  III. 


The  heading  of  the  projectile,  h,  is  in  the  direction  of  its  air  velocity.  It  is  also  de 
fined  in  Chapterlll: 


I  | 


'  •* 

In  terms  of  h  and  D,  the  aspect  angle  de,„  of  the  projectile  as  seen  from  the  radar  is 
e„p.cos-\b.h)  (4.21) 


Note  that  0  <  6ttp  <  n.  For  simplicity,  the  projectile  is  considered  as  a  circular  cone  on 
top  of  a  circular  cylinder,  as  shown  in  Figure  24  on  page  37.  The  height  of  the  cylinder 
is  414.4  mm.  The  height  of  the  cone  is  523.6  mm.  The  diameter  of  the  cylinder  and  the 
base  of  the  cone  is  155  mm. 
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Figure  24.  Approximate  Shape  of  155mm  Projectile. 


The  vertex  angle  of  the  cone,  2a,  is  given  bv: 


a  =  tan 


2  x  523.6 


(4.22) 


=  0. 1469  rad. 


(4.23) 


=  8.42° 


(4.24) 


The  physical  cross-section  of  the  cone  is  the  projected  area  of  the  projectile  along  D.  the 
radar  pointing  direction,  onto  a  plane  perpendicular  to  D.  This  cross-sectional  area  will 
be  discussed  below  according  to  whether  |  D  •  h  |  <  cos  a  or  |  D  •  h  |  >  cos  a. 

1.  Case  1 

As  |  D  •  h  |  ^  cos  a,  the  projection  of  the  projectile  along  D  is  shown  in 
Figure  25  on  page  3S. 
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Figure  25.  Radar  Cross  Section  Case  I. 


The  total  RCS  value  is 

a ,  =  na 2  \  D  •  h  |  +  lab  |  D  x  h  \  (4.25) 

where  a  in  the  radius  of  the  cylinder  and  b  is  the  height  of  the  cylinder. 

2.  Case  2 

As  1  D  •  /;  |  <  coso:.  the  total  RCS  is  given  by  equation  (4.20)  (Figure  26  on 
page  39). 


Figure  26.  Radar  Cross  Section  Case  2. 


Ol  =  <7)  +  a  2  +  ct3  +  <74 


M26) 


ZA  tan  a(  tan'fl  4-  tan^a V ' 2(  tan20  -  tan‘a)w(P  •  A") 

tan‘0 


(4.27) 


ol.  , ,  _  C0S  ®a'P  s3/2 

2  v  *  2  * 

sin  6asp  cos  a 


(•4.28) 


2,  _■  .-i,  tan  a  ,  tan  a 
°  ‘  |  tan  6asp  |  |  tan  oasf_ 


1 - taalg_  |  |  f)./J 

tan  0asF 


(4.29) 


cr3  =  2fl6  !  D  x  h 


(4.30) 


A  A 

2  I  D  .  h  I 
er4  =  na  - - - 


,V  = 


D  cos  a  —  h  cos  6 


asp 


x  sin^0av  cos4a  +  sinJa  cos ~0a!p 


(4.31) 


(4.32) 


where  L  is  the  height  of  the  cone  and  .V  is  the  normal  vector  of  the  surface  Ar  The  de¬ 
tailed  derivation  of  the  RCS  equations  is  given  in  Appendix  B. 

From  the  above  results,  the  RCS  value  of  the  projectile  in-flight  can  be  pre¬ 
dicted.  It  is  shown  in  Figure  27  and  Figure  28  on  page  41.  From  above  figures,  it  can 
be  seen  that  the  maximum  RCS  is  about  0.09  m2  and  the  minimum  RCS  is  about 
0.02  m2.  Since  the  aspect  angle  of  the  projectile  changes  very  fast  right  after  it  is  fired, 
the  RCS  value  varies  rapidly  at  the  beginning  of  the  trajectory. 


Figure  27.  The  In-Flight  Projectile  Radar  Cross  Section  vs.  Aspect  Angle. 
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D.  PROPAGATION  EFFECTS 

Propagation  of  radio  waves  through  the  Earth's  atmosphere  is  significantly  affected 
by  the  propagation  medium.  The  presence  of  the  Earth  complicates  the  situation  in 
three  ways:  (a)  the  surface  of  Earth  scatters  radiation  and  produces  an  interference  pat¬ 
tern,  (b)  the  Earth  casts  a  shadow  and  gives  rise  to  diffraction  phenomena,  and  (c)  the 
Earth's  atmosphere  is  inhomogeneous  and  produces  refraction  effects  [Ref.  13].  It  is 
customary  to  group  all  these  factors  into  a  single  quantity  called  the  propagation  factor 
(F).  In  this  thesis,  the  Engineer's  Refractive  Effects  Prediction  system  (EREPSj  [Ref. 
14]  is  utilized  to  compute  the  propagation  factor  based  on  local  environment. 

EREPS  uses  ray  tracing  for  radio  wave  propagation  within  the  horizon.  It  includes 
the  bending  effect  of  the  refractive  index  variation  of  air  and  the  divergence  factor  when 
the  ray  is  reflected  from  the  Earth.  The  field  at  the  target  is  the  sum  of  the  direct  ray 
and  the  reflected  ray,  resulting  in  interference  peaks  and  nulls  which  become  significant 
when  the  target  falls  near  the  horizon  (Figure  29). 


/TARGET 


Figure  29.  Reflection  and  Interference  from  a  Spherical  Earth. 
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For  the  over-the-horizon  region,  the  waveguide  mode  approach  is  adopted  by 
EREPS.  The  Earth's  atmosphere  is  modeled  as  a  leaky  waveguide.  The  Earth’s  surface 
is  the  lower  boundary  of  the  guide.  VI  odes  of  high  attenuation  along  the  horizontal  di¬ 
rection  are  called  leaky  modes.  Modes  of  low  attenuation  are  called  trapped  modes.  In 
x-band,  the  effect  of  the  evaporation  duct  is  significant.  Trapped  modes  occur  when  the 
evaporation  duct  height  is  large  or  at  least  comparable  to  the  radar  wavelength  (see 
Figure  30  on  page  44)  (Ref.  15]. 

1.  Computation 

A  typical  EREPS  output  is  given  in  Figure  30  on  page  44.  The  local  annual 
average  meteorological  conditions  for  propagation  are  available  [Ref  13].  The  following 


parameters  are  used  : 

1.  Evaporation  duct  height  :  13.1m 

2.  k  (effective  Earth  radius  factor)  :  1.45 

3.  Wind  speed  :6ms 

4.  Absolute  humidity  :  lOg  m: 

5.  Surface  refractivity  :  370. 


Note  that  the  projectile  is  at  a  particular  height  when  it  arrives  at  a  range.  This 
computation  has  to  be  repeated  for  eacli  point  along  a  projectile  trajectory.  This  is  done 
for  several  trajectories  at  various  ranges  (Figure  32  on  page  45  and  Appendix  C  ).  Also 
note  that,  in  Figure  32  on  page  45.  the  ground  range  is  given  a  negative  sign  because  it 
is  measured  from  the  point  of  impact  toward  the  radar. 

2.  Summary  of  Results 

a.  Signal  Strength  Enhancement 

From  Figure  32  on  page  45,  it  can  be  seen  that,  with  ducting,  which  is 
present  most  of  the  time  at  the  test  site,  the  return  signal  will  be  substantially  enhanced 
from  the  corresponding  free  space  value.  (Note  that  the  computation  is  carried  out  for 
a  target  at  1  m  height,  i.e..  when  the  projectile  is  near  its  point  of  impact.  The  applicable 
range  for  the  projectile  of  interest  lies  between  27  to  28  km.)  Thus  ducting  is  helpful  in 
enhancing  the  quality  of  the  recorded  full  range  tracking  data. 

b.  Strong  Interference 

As  shown  in  Figure  32  on  page  45.  the  return  signal  will  exhibit  strong  in¬ 
terference  when  the  projectile  reaches  about  250  m  within  the  point  of  impact.  The  more 
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than  ±15  dB  fluctuation  makes  it  difficult  to  track  the  projectile  in  this  region  using  the 
simple  CW  system  proposed  in  this  thesis.  On  the  other  hand,  since  this  strongly  fluc¬ 
tuating  region  is  within  a  few  hundred  meters  from  the  impact  point,  this  small  portion 
of  the  trajectory  can  be  estimated  with  good  accuracy.  Hence  it  is  desirable  to  select  a 
'stop-processing  point'  at  which  the  CW  real-time  radial  velocity  and  range  computa¬ 
tions  will  be  halted.  The  radar  antenna  will  be  guided  by  the  predicted  trajectory,  and 
the  return  video  signal  will  still  be  recorded. 


Figure  30.  Waveguide  Mode  with  Ducting  Height  and  Antenna  Height. 
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Propagation  factors  (dB) 
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re  32.  Propagation  Factors  vs  Ground  Range  from  Impact  Point. 


3.  Stop-Processing  Point 

It  should  be  noted  that  EREPS  does  not  take  into  account  the  antenna  radi¬ 
ation  pattern.  Since  the  antenna  to  be  used  have  a  beamwidth  of  1°  only,  interference 
between  the  direct  ray  and  the  reflected  ray  will  not  be  important  until  the  antenna  ele¬ 
vation  angle  is  lowered  to  within  0.5°  above  the  horizontal  direction  (Figure  33  on  page 
46). 


Figure  33.  Narrow  Beamwidth  Effect  on  the  Spherical  Earth  Surface. 


This  criterion  can  be  used  to  determine  the  stop-processing  point  as  follows:  when  the 
system  records  an  antenna  elevation  angle  falling  below  0.5C,  the  square-law  detector 
and  the  comparator  will  be  switched  off.  The  servo  will  steer  the  antenna  according  to 
the  predicted  trajectory  smoothly  toward  the  projectile  point  of  impact.  Characteristics 
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for  two  such  stop-processing  points  for  a  radar-to-gun  distance  of  500  m  and  two  dif¬ 
ferent  radar  heights  of  10  m  and  30  m  respectively  are  displayed  in  Table  1  on  page  47. 
Note  that  these  point  are  close  to  where  the  propagation  factor  starts  to  fluctuate 
strongly  as  predicted  by  EREPS. 

Since  the  processing  of  real-time  data  is  stopped  before  contributions  from  the 
reflected  ray  become  important;  in  what  follows,  the  propagation  factor  will  be  consid¬ 
ered  as  unity. 


Table  1.  STOP  PROCESSING  POINT. 


Antenna  height  (m) 

10 

30 

Radar  to  gun  distanced rm 

500 

500 

Projectile  height  (mi 

30S 

305 

Slant  range  (kmt 

27.53 

27.53 

Impact  range  (km) 

27.69 

27.57 

Ground  range  to  impact  point  (m) 

15S 

157 

Elevation  angle  (degree) 

0.49 

0.45 

Elevation  angle  at  impact  point  (degree) 

-0.15 

0. 1 9 

Elevation  angle  difference  (degree) 

-0.64 

-0.64 

Azimuth  angle  (degree) 

36.5 

36.5 

Azimuth  angle  at  impact  point) degree) 

36.9 

36.9 

Azimuth  angle  dillcrencc  to  impact  point(dcgree) 

0.39 

0.39 

In-flight  time  (sec) 

S6.4S 

S6.4S 

E.  POWER  REQUIREMENT 

From  equation  (4.1).  the  power  requirement  for  tracking  up  to  the  stop-processing 
point  is  shown  in  Figure  34.  In  order  to  reduce  the  required  transmitter  power  re¬ 
quirement,  a  multiple  filter  design  is  proposed  for  this  system. 

The  multiple  filter  design  is  based  on  the  characteristics  that,  after  an  initial  rise,  the 
doppler  frequency  decreases  with  increasing  slant  range  (Figure  35  on  page  49).  By 
switching  to  video  filters  of  narrower  bandwidths  at  longer  ranges,  the  coherent  inte¬ 
gration  time  of  the  radar  is  increased  hence  the  power  requirement  can  be  lowered. 
Combining  the  fixed  range  independent  parameters  in  the  radar  range  equation,  the 
power  requirement  can  be  expressed  as  : 


P,  =  (3.7339  x  10_2i)  x  ( 


(4.33) 


=  C  x  (  —  )maxfd„ 


(4.34) 


C  =  4.4807  x  10" 


where  Z?,  is  1.2/". 


Figure  34.  Power  Requirement  for  Full  Range  Tracking  without  Multiple  Filters. 


Since  fd.a  and  R  are  functions  of  the  trajectory  (Figure  35  on  page  49.  and  Figure  36 
on  page  50).  if  the  tracking  range  is  divided  into  three  regions  1,  2,  and  3,  and  one  filter 
is  used  for  each  region,  the  transmitter  power  requirement  becomes  : 


—  C(  ~a~  ): .max  {fjh.max 


~  C(  G  )3 tmax  if 3, max 

where  the  constant  C  is  the  same  as  that  of  equation  (4.34). 


(4.36) 


(4.37) 


(4. 38) 


Figure  35.  Doppler  Frequency  vs.  Slant  Range. 


The  power  requirement  written  in  terms  of  that  of  the  single  filter  case  is: 


P,- 


R 4 

^  n  )maxfdmax 

c„ 


(4.39) 


where  C..  the  power  reduction  factor,  is  the  factor  by  which  the  power  requirement  is 
reduced. 

In  order  to  find  the  optimal  switch-over  points  for  the  filters,  the  logarithm  of/, 
versus  7?4/<7,  normalized  to  (/)m„  and  (/?J/a)m„  respectively,  is  plotted  (Figure  3r  on  page 


Assume  that  the  first  switch-over  point  is  located  at  and  the  second  one  is 

located  at  (jr^jj).  These  values  can  be  determined  from  Figure  37  on  page  51  by  ob¬ 
serving  that,  at  the  optimal  switch  points  : 


II 

+ 

(4.40) 

~  >2 

(4.41) 
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The  optimal  locations  so  obtained  are  listed  in  Table  2  on  page  52.  The  power  re¬ 
duction  factor  is  CP  =  100621!  =  4.1S6.  as  seen  from  Figure  37  on  page  51. 

The  power  requirement  as  a  function  of  slant  range  when  three  filters  are  used  is 
plotted  in  Figure  38  on  page  52.  The  antenna  height  is  10  m  and  the  radar-to-gun  dis¬ 
tance  is  500  m.  The  required  power  is  37  Watts.  A  conventional  CW  tracking  system 
will  track  to  the  summit  point  only.  This  point  is  located  at  a  range  of  region  I  and  re¬ 
quires  a  transmitter  power  of  19.5  Watts.  Almost  full-range  tracking  is  possible  by 
doubling  the  conventional  power  through  the  use  of  multiple  filters. 


Figure  37.  Log(Normalized /<)  vs.  Log( Normalized  R*lo). 
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Figure  38.  Power  Requirement  for  Trajectory  Tracking  with  Multiple  Filters. 


Table  2.  MULTIPLE  FILTERS  SWITCH  POINTS 


Parameters 

region  1 

region  2 

region  3 

stop  track¬ 
ing  point 

Slant  range  (km) 

16.06 

27.53 

27.53 

Maximum/.  (kHz) 

22.03 

12.72 

... 

Video  bandwidth 
(kHz) 

63.90 

26.44 

15.27 

... 

Radar  cross  section 
( m: ) 

0.03S5 

0.04  3  S 

0.0903 

Required  power 
(Watts) 

36.29 

36.23 

36.29 

... 
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V.  CONCLUSIONS 

It  is  demonstrated  in  this  thesis  that  full-range  projectile  tracking,  with  real-time 
track  processing  stopped  at  a  range  within  3  60  m  of  the  projectile  point  of  impact,  can 
be  achieved  if  the  transmitter  output  power  is  37  watts  or  higher  and  a  set  of  three 
switchable  filters  are  used.  If  a  single  filter  is  used,  an  output  power  of  160  watts  is  re¬ 
quired.  Either  configuration  can  be  implemented  with  a  commercially  available,  small 
size,  air-cooled  RF  amplifier  which  can  be  mounted  near  the  transmitting  antenna. 

The  tracking  relies  on  the  predicated  trajectory.  Although  real-time  doppler  infor¬ 
mation  is  available  from  the  radar,  it  is  used  only  to  confirm  that  the  projectile  is  on 
track.  A  algorithm  to  utilize  this  real-time  information  for  updating  the  predicted  tra¬ 
jectory  is  desired.  Kalman  filtering  appears  to  be  a  promising  technique  which  deserves 
further  investigation. 


APPENDIX  A.  PROJECTILE  TRAJECTORY  SIMULATION  AND 

RANGE  ANALYSIS 


clear 
t2=clock; 
format  long  e; 

X  The  first  capital  letter  designations  are:  P  for  the  projectile; 

X  G  for  the  gun  (muzzle);  R  for  the  radar. 

X  The  numerals  1,2,3  are  for  the  radial,  theta  and  phi  components  in 
X  an  inertial  spherical  coordinate  system;  4  is  for  the  radial  component 
X  in  the  cylindrical  coordinate  system  associated  with  the  spherical 
X  coordinate.  The  origin  is  at  the  center  of  the  Earth;  the  polar  (z) 

X  axis  is  the  axis  of  spin  of  the  Earth. 

X  East  and  South  designations  are  local  quantities  in  a  frame  rotating 
X  with  the  Earth 


X  constant  declaration 
EarthT=86164. 06; 
EAV=2*pi/EarthT; 
a=( 63569 12+6378533 )/2; 
g=9.  81; 


%  sidereal  period  of  Earth  spin 
X  Earth  angular  velocity 
X  approximate  Earth  radius 
%  gravitational  acceleration 


RLATI=24+36/ 60+38/ 3600; 

RL0NGI=12 1+5 2/60+27/ 3600; 
f=10525£6; 
lambda=3E8/f ; 

?0input  parameters 
! de 1  ww 
diary  ww 

RB=input( ' radar  to  gun  base  distance(m)  =?  ') 

X  RBaz=input( '  radar  to  gun  direction  from  east=?  ' )*pi/180 
RBaz=0; 

X  MN=input( ' Measurement  base  number  setting  =?  ') 

MN=1000; 

Rh=input( ' Radar  height(m)=?  ') 

X  Gel=input( ' firing  elevation  angle( degree)=?  ' )/180*pi 
Gel=45/180*pi; 

X  Gaz=input( 1  firing  azimuth  angle(degree)=?  ')/18G,vpi 


X  latitude  24o36'38"  (radar  site) 

X  longitude  121o52'27"  (radar  site) 
X  radar  transmitted  frequency  Hz 
X  radar  transmitted  wavelength  m 


Gaz=45/180*pi; 

diary 

answer=input( ' If  you  want  reenter  parameter  please  enter  1  else  enter 
if  answer=l 
return 
xxa 


end 

B 1=5 ;  X  gun  barrel  length 


A=pi*( 0, 155/2 )**2; 
m=45. 4; 

FAC=0. 5*A/m; 


X  projectile  cross  section  area 
X  projectile  mass  ( kg) 

X  coefficient  for  computing  drag 
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MB=MN*( lambda/ 2 ) ; 


%  measurement  base 


%  Radar  Location  at  t=0 
Rl=a+Rh; 

R2=( 90-RLATI )*pi/ 180; 

R3=0; 

R2s=sin(R2); 

R2c=cos(R2); 

R4=Rl*R2s; 

Rz=Rl*R2c; 

RV=R4*EAV; 

%  Firing  Direction  and  Speed 
PVo=850; 

Gels=sin(Gel ); 

Gelc=cos(Gel); 

Gazs=sin( Gaz) ; 

Gazc=cos( Gaz); 

%  Gun  (muzzle)  Location  at  t=0 
Gh=Bl*Gels; 

Gl=a+Gh; 

Gl=Bl*Gelc; 

RGs=-RB"'sin(  RBaz)  -Gl*Gazs; 

RGe=RB''rcos(RBaz)+Gl*Gazc; 

G2=R2+RGs/Gl; 

G4=Gl*sin(G2); 

G3=R3+RGe/G4; 

Gz=Gl*cos(G2); 


%  the  inertial  z-x  plane 

%  z-coordinate  of  radar 
%  radar  velocity 

%  muzzle  speed  of  projectile 

%  muzzle  height 

%  distance,  south  of  radar 
%  distance,  east  of  radar 


«  Displacement  Transformation 
G2c=cos( G2 ) ; 

G2s=s in( G2 ) ; 

DTG=( 0 ,GAZC ,GAZS; 0 , -GAZS ,GAZC; 1 ,0 , 0-*-G2C , 0 ,G2S; 0 , 1 , 0; -G2S,0,G2C); 

%  Projectile  Initial  Location 
P1=G1; 

Pln=Gl; 

T’2n=G2; 

P3n=G3; 

P4n=G4; 

%  Projectile  Initial  Velocity 
PVln=PVo*Gels; 

PVg=PVo*Gelc; 

PVe=PVg*Gazc; 

PA V  2n= - ( PVg*Gazs ) / G 1 ; 

PAV3n=PV'e/G4+EAV; 


«  velocity  along  ground 
%  velocity,  east 
%  angular  velocity,  theta 
%  angular  velocity,  phi 


%  Drag  Force  Considerations 

%  Wind  Profile:  assumed  constant  in  this  program. 

%  The  geostrophic  wind  should  be  measured  and  the  Ekman  spiral 
model  can  be  adopted  for  a  better  approximation  if  the  true 


%  wind  profile  is  not  available. 

%  Wind=input( ' Enter  the  wind  speed  m/sec  at  muzzle  height') 

Wind=6;  %  wind  velovity  at  muzzle  height 

%  WindBearing=input( ' Enter  the  true  wind  bearing  in  degrees') 
WindBearing=135;  %  clockwise  from  north 

Waz=( (90-WindBearing)+180)*pi/ 180; 

We=Wind*cos(Waz);  %  wind  velocity,  east 

Ws=-Wind*sin(Waz);  %  wind  velocity,  south 

%  Temperature  Profile  (assumed  linear) 

%  Th=input( ' Enter  the  height  (m)  at  which  temperature  is  measured') 
Th=10; 

%  TREF=INPUT( 'ENTER  THE  MESURED  TEMPERATURE  IN  CELSIUS') 

Tref=15; 

TrefK=273. 16+Tref; 

Thl=Th+a; 

S  =-0. 0065; 

%  TP=TREFK+S* ( P 1 -TH 1 ) ; 

%  AIR  DENSITY 

%  DH=INPUT( 'ENTER  THE  HEIGHT  (M)  AT  WHICH  AIR  DENSITY  IS  MEASURED') 
DH=GH; 

%  DA0=INPUT( 'ENTER  THE  MESURED  AIR  DENSITY  IN  KG/CUBIC  METER') 

Dao=l. 225;  %  air  density  at  sea  level 

TDao=TrefK+s*( Dh-Th) ; 

ss=287. 05;  %  gas  constant  N-m/kg-K 

KC^l. 4;  %  Cp/Cv  for  air 

RAg=-l-g/(s*ss);  %  calculation  factor 

%  DA=DAO*(TP/TDAO)**RAG;  %  AIR  DENSITY  AT  PROJECTILE  HEIGHT 

%  Speed  of  Sound 
KR=KC*ss; 

%  Vah=sqrt(KR‘vTP);  %  speed  of  sound  at  PI 

%  curve  fitting  for  Drag  coefficient  experimental  data  input) 
x=( 0,0. 05,0. 15,0. 2,0. 23,0. 25,0. 3,0. 35 ,0. 4,0. 45,0. 5,0. 55,0. 6,.  .  . 

0.  65,0.  7,0.  75,0.  825,0.  875,0.  98,1.  0,1.  05,1.  1,1.  2,1.  3, 1.4,1.  5,.  .  . 

1.  75,2,2. 5,2. 9); 

y=(0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,0. 16,.  .  . 

0. 16,0. 16,0.  16,0.  16,0.  162,0. 17,0. 205,0. 225,0. 27,0. 282,0. 29,. . . 

0. 2875,0. 285,0. 275,0. 25,0. 23,0. 195,0. 17); 

cds=spline(x,y);  %  spline  curve  fitting  function 

n  get  spline  coefficient 

%  Trajectory  Computation  Begins 

t(l)=0;  %  initial  time 

n=0; 

%*vr***iv********tfrfr****  main  loop  begin  ******************************* 
while  Pln>a,  n=n+l; 

%  Update 
P l=Pln; 

P2=P2n; 
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%  TEMPERATURE  LAPSE  RATE  K/M 
%  AIR  TEMPERATURE  AT  THE  PROJECTILE 


P3=P3n; 

P4=P4n; 

PVl=PVln; 

PAV2=PAV2n; 

PAV3=PAV3n; 

PV2=PAV2*P1; 

PV3=PAV3*P4; 

PVe=( PAV3-EAV)*P4; 

Height(n)=Pl-a; 

%  velocity  transformation  matrix 

PR3=P3-EAV*t(n);  %  location  relative  to  radar 

PG3=PR3-G3;  %  location  relative  to  gun 

P2s=sin(P2); 

P2c=cos(P2); 

PR3s=s in( PR3) ; 

PR3c=cos( PR3); 

PG3s=sin(PG3); 

PG3c=cos(PG3); 

VTR1 l=P2s*PR3c; 

VTR21=P2s*PR3s; 

VTR31=P2c; 

VTR=(VTRll,F2c*PR3c,-PR3s; VTR21,P2c*PR3s,PR3e; VTR31,-P2s,0); 

VTG1 l=P2s*PG3c; 

VTG21=P2s,vPG3s; 

VTG31=P2c; 


%  Relative  Location 


PRXYZ=( P 1* - VTR 1 1 ;  VTR21;  VTR31)-(R4;  0;  RZ1; 

S LR=N ORM ( PRXY Z ) ;  %  SLANT  RANGE  FROM  RADAR 

SRR(n)=SLR; 

®_  '.V  Vr  Vr  Vr  Vr  Vr  Vr  V:  Vr  V :  Vr  *  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  W  Vr  Vr  Vr  Vr  Vr  Vr  -.V  V:  Vr  *  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  '.V  Vr  Vr  Vr  Vr  Vr  -.V  *  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr 


DR( :  ,n)=PRxvz/SLR;  ?i 

Dz=DR ( 1 , n ) *R2  s+DR ( 3 , n )*R2  c ;  % 

Dn=DR( 3 ,n)*R2s-DR( 1 ,n)*R2c; 
De=DR(2,n);  % 

traz(n)=atan2(Dn ,De)* 180/pi;  % 

trel(n)=asin(Dz)*180/pi;  % 

^Vr  Vr  VrVr  Vr  Vr Vr  V?  Vr  Vr VrVr  V?  Vr  Vr  Vr  VrVr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  VrVr  VrVrVr 


direction  relative  to  radar 
tracking  direction,  elevation 
tracking  direction,  north  comp, 
tracking  direction,  east  component 
azimuthal  angle 
elevation  angle 

‘Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  *  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vt 


PGXYZ=DTG*(P1*-VTG11;  VTG21;  VTG31)-(G4;  0;  GZ) ) ; 
SL=norm( PGxyz) ;  %  slant  range 

SRange( n)=SL; 

GRange( n)=a*acos ( ( P 1**2+G 1**2 - SL  2 ) / ( 2*Pl*G 1 ) ) ; 

%  ground  range 


PGx=PGxyz( 1); 

PGz=PGxyz( 3); 

FRange(n)=a*atan(PGx/(Gl+PGz) );  %  range  in  firing  plane 

CRange(n)=PGxyz( 2);  %  cross  range  from  firing  plane 


^Relative  Velocitv 
PV=( PV1;  PV2;  PV3); 

PVR=VTR*PV-(0;  RV;  0); 

VR( :  ,n)=PVR;  %  velocity  relative  to  radar 

VRr(n)=PVl; 


Time  Interval  and  Doppler  effect  Estimation 


%  Doppler  velocity 


VDopp=PVR ' *PRxy z/ SLR; 
Vdopl(n)=VDopp; 
fdl(n)=2,vVDopp/ lambda; 
dtn=MB/VDopp; 
t(n+l)=t(n)+dtn; 

%  Frictional  Acceleration 

TP=TrefK+s*( PI -Thl ) ; 

Da=Dao*( TP/TDao ) **RAg; 
Vah=sqrt(KR*TP); 

PW=(  PV1;  PV2-WS;  PVE-WE); 
PVWM=norm(PW); 

PWMM=PVWM/Vah; 

PH=PVW/PWM; 

HeadingR( : ,n)=VTR*PH; 

CD=ppva  1  ( cds  ,  PVWMM ) ; 

FA=  -  F  AC*D  a*CD*P  WM*PVW ; 

%  Location  and  Velocity  Prediction 


%  Doppler  Frequency 
%  square-law  doppler  detector 


%  air  temperature  at  the  projectile 
%  air  density  at  projectile  height 
%  speed  of  sound  at  PI 
%  PROJECTILE  AIR  VELOCITY 


%  projectile  heading 
°o  heading  as  seen  from  radar 

%  drag  acceleration 


PA1=( PV2**2+PV3**2 ) /P1+FA( 1 ) -g*( a/Pl )  2; 

PAA2=P2s*P2c*PAV3**2+FA( 2) /PI; 

PAA3=FA( 3)/P4; 

PVln=PAl*dtn+PVl; 

Pln=0.  5 * ( P V 1 +PV 1  n ) ,vd tn+P  1 ; 
PAV2n=(PAA2*dtn+PAV2)*(Pl/Pln)**2; 

P2n=0. 5*( PAV2+PAV2n)*dtn+P2; 

P4n=Pln*sin(P2n); 

PAV3n=(  PAA3‘-Vdt n+P AV3 ) *(  P4 / P4n )  ** 2 ; 

P3n=0. 5*( PAV3+PAV3n)*dtn+P3; 
end 

a*************  end  of  main  loop  *************** 
%  Estimate  Point  of  Impact 
Ih=Pl -a; 

IA=-PA1; 

IV=-PV1; 

dt=(sqrt(IV**2+2*Ih*IA)-IV)/IA; 

N=n+1; 

t(N)=t(n)+dt; 

Pln=a; 

PVln=PAl*dt+PVl; 

PAV2n=(PAA2*dt+PAV2)*(Pl/Pln)**2; 

P2n=0. 5*( PAV2+PAV2n)*dt+P2; 

P4n=P ln*s in( P2n ) ; 

PAV3n=(  PAA3‘-''dt+PAV3)*(  P4/P4n)**2; 

P3N=0. 5*(PA\3+PAV3N)*DT+P3; 


%  Update 
Pl=Pln; 
P2=P2n; 
P3=P3n; 
P4=P4n; 
PVl=PVln; 
PAV2=PAV2n; 
PAV3=PAV3n; 
PV2=PAV2*P1; 
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PV3=PAV3*P4; 

Height(N)=0; 

%  velocity  transformation  matrix 

PR3=P3-EAV*t(N);  %  location  relative  to  radar 

PG3= PR3-G3;  %  location  relative  to  gun 

P2s=sin(P2); 

P2c=cos( P2 ) ; 

PR3s=sin( PR3 ) ; 

PR3c=cos( PR3 ) ; 

PG3s=sin(PG3); 

PG3c=cos(PG3); 

VTRll=P2s*PR3c; 

VTR21=P2s*PR3s; 

VTR31=P2c; 

VTR=( VTR 1 1 , P2C*PR3C , -PR3S;  VTR2 1 , P2C*PR3S , PR3C;  VTR3 1 , -P2S , 0 ) ; 

VTGll=P2s*PG3c; 

VTG21=P2s*PG3s; 

VTG31=P2c; 


%  Relative  Location 

PRXYZ=( P1*-VTR11;  VTR21;  VTR311 -( -R4;  0;  RZ)-; 

SRR(  N')=norm(  PRxyz); 

*********************  Radar  parameters  registrat ion  ****************** 


DR(: , N ) =PRxy z / SRR ( N ) ; 

Dz =DR  ( 1 ,  N  ) *R  2  s  +DR  ( 3 ,  N  )  ■ *R  2  c ; 
DN=DR( 3 ,N)*R2S-DR( 1 ,N)*R2C; 
De=DR( 2 ,N) ; 

traz(N)=atan2(Dn,De)* 180/pi; 
t re l(N)=asin(Dz)* 180/pi; 

*  PGxyz=DTG*(Pl*[ VTG11;  VTG21; 
SL=norm( PGxyz ) ; 

SRange(N)=SL; 

GRange( N )=a*acos ( ( P 1**2+G 1**2 -SL 


location  relative  to  radar 
tracking  direction,  elevation 
TRACKING  DIRECTION,  NORTH  COMP, 
tracking  direction,  east  component 
azimuthal  angle 
elevation  angle 

W  V,  V-  J#  Jb  Jb  bV  Jb  bib  ■'* b'b  Jb  Jb  Jb b'b  b'b 

VTG31]  -[  G4;  0;  Gz]  ); 

%  slant  range 


PGx=PGxyz( 1); 

PGz=PGxyz( 3 ) ; 

FRange(N’)=a*atan(PGx/(Gl+PGz)); 
CRange(N)=PGxyz(  2); 


2)  /  ( 2*P1*G1 ) ) ; 

%  ground  range 


%  ground  range  in  firing  plane 
%  cross  range  from  firing  plane 


%  Relative  Velocity 
PV=[ PV1;  PV2;  PV3] ; 

PVR=VTR*PV-(  0;  RV;  0]; 

VDopp=PVR ' *PRxyz/SLR;  %  Doppler  velocity 

Vdop 1 ( N ) =VDopp; 

fdl(Nj=2*VDopp/lambda;  %  Doppler  Frequency 

VR( : ,N)=PVR; 

VRr( N)=PV1; 

PVW=( PV1;  PV2-WS;  PVE-WE); 

HeadingR( :  ,N)=VTR*(  PW/norm(  PW) ) ; 
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for  i=l:N; 

VRM( i)=norm(VR(: ,i)); 

VRrf( i)=VRr( i)/VRM( i) ; 

VR t ( i ) =( 1 - VRr  f ( i ) ) *VRM( i ) ; 
beamD( i)=l*pi/180*SRR( i); 
end 

EXETIMEl=etime( clock, t2) 

J**fnV**Vr***)W:*Sr**  OUt  put  data  **************************** 

diary  ww 
EXETIME1 
MN=MN 
ng=N’ 

MB=MB 

tg=t(N) 

r=Pl 

dlatitude=( P2-G2)* 180/pi 
dlongitude=PG3* 180/pi 
Rg=GRange(N) 

FRg=FRange(N) 

CRg=CRange(N) 

CRgmax=max( CRange ) 
maxHeight=max( Height) 
maxfd=max( fdl ) 
maxrange=SRR(N) 

VDOPmin=min(  Vdopl ) 

VDOPmax=max( Vdopl ) 
diary 

********************  SAVE  AND  LOAD  DATA  FOR  NEXT  CALCULATION'***** 
50010. MAT 

50010  MAXRANGE  EXETIME1  DR  HEADINGR  SRR  HEIGHT  T  N  FD1  MB 
X50010. MAT 

X50010  TRAZ  TREL  VD0P1  VRM  VRRF  VRT  VR  VRR 

Execute  other  program 


IV  -'-.u-u, 
A>" 

!  DEL 
SAVE 
!  DEL 
SAVE 


%PR 

"oPRTRAK 
%XX0 
%RCSA 
GRANGE 3 
%TDATA 


v******v? 


-/o.i  ..  ..  ..  ...r-.,. ..  ..  ..  ..  ..  ..  .<  ..  ..  ..  .. ---output 

! del  GH. met 

plot(GRange .Height) ,grid; 
title( 'PROJECTILE  TRAJECTORY'); 
xlabel( ' ground  range  (m)'); 
ylabelf 'projectile  height  (m)'); 
meta  GH 
«prtsc 
%pause 

! del  GC. met 

plot(GRange, CRange) ,grid; 
tit le( 'PROJECTILE  TRAJECTORY’); 
xlabel( ' ground  range  ( m ) ' ) ; 
ylabel( ' Cross  range  ( m ) ' ) ; 


SRR  BEAMD 

******************** 
PRINT  OUT  TRAJECTORY 
TRACKING  PARAMETERS  &  OUTPUT 
RCS  &  RANGE  ANALYSIS  PROGRAM 
RCS  SUBFUNCTION 
POWER  REQUIREMENT 
PROPAGATION  FACTOR  CURVE  FITTING 

of  traject 0 r y********************* 
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meta  GC 
%  pause 
»  prtsc 

! del  RH. met 

p lot ( SRR , Height , 'o' ) , grid; 
t it le( 'PROJECTILE  TRAJECTORY'); 
xlabel( ' Slant  Range  (m)'); 
ylabel( 'projectile  height  (m)'); 
meta  RH 
%  prtsc 
pause 

*>  ! del  TH.  met 
plot(t, Height) .grid; 
xlabel( ' flight  time  (sec)'); 
ylabel( ' projectile  height  (m)'); 
meta  TH 
%  pause 
?oprtsc 

! del  TSGFC.  met 

plot(t,SRange,  t.GRange,  t, FRange,  t .CRange) .grid; 
xlabel( ' f 1 ight  time  (sec)'); 

ylabel( ' slant ,  ground,  firing  plane  &  cross  ranges  (m)'); 
meta  TSGFC 
/ipause 
prtsc 

! del  SGFCT.  met 

plot(SRange,t,  GRange,t,  FRange, t,  CRange ,t) ,grid; 
xlabel( 'slant ,  ground,  firing  plane  &  cross  ranges  (m)'); 
ylabel( ' flight  time  (sec)'); 
meta  SGFCT 
%pause 
?oprtsc 

! del  TYR.  met 
plot( t , VRr , ' +' ) ,grid; 
xlabel( ' flight  time  (sec)'); 
ylabel( ' Vertical  Speed'); 
meta  TVR 
%pause 
%prtsc 
! del  TVH.  met 
plotf  t , VRt , '+' ) .grid; 
xlabel( 'Flight  time  (sec)'); 
ylabel( ' Horizontal  Speed'); 
meta  TVH 
%pause 
prtsc 

! del  FH. met 

p lot (FRange .Height) ,grid; 

xlabel( ' ground  range  in  firing  plane  (m)'); 
ylabel( 'projectile  height  (m)T); 
meta  FH 
%pause 
%prtsc 
! del  FC. met 

plot (FRange .CRange) ,grid; 
xlabel( 'Firing  plane  Range  (m)'); 
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ylabel( ' Cross  range  (m)'); 
meta  FC 
%pause 
prtsc 

!del  GV.  met 
plot(GRange ,VRM) .grid; 
xlabel( ' Ground  Range  (m)'); 
ylabel( ' Projectile  Speed  (m/sec)'); 
meta  GV 
%pause 
%prtsc 
! del  GRV. met 
plot(GRange ,VRrf ) .grid; 
xlabel( ' Ground  range  (m)'); 

ylabel( 'Ratio  of  radial  velocity  to  speed'); 
meta  GRV 
%pause 
%prtsc 

*****************  TRAKING  parameters  calculation  and  output  ****** 

%****************  Tracking  angular  velocity  ************ 
avell=diff(trel). /diff(t); 
avazl=diff(traz). /diff(t); 

********************  PLOT  ****************************** 
for  n=l : N-l 
ta(n)=t(n) ; 
end 

!del  AVAEt.  met 
plot(ta  .avell , '  +  ' ) .grid; 

title( 'RADAR  TRACKING  ELEVATION  ANGULAR  SPEED  '); 
xlabel( ' In-flight  time  (sec)'); 

ylabel( ' Elevation  Angular  velocity  (degree/sec)'); 
meta  AVAEt 
%  pause 
prtsc 

! del  AVAZt.  met 

plot( ta , avazl , ' +' ) .grid; 

title( 'RADAR  TRACKING  AZIMUTHAL  ANGULAR  SPEED’); 
xlabel( ' In-flight  time  (sec)'); 
ylabel( ' Azimuthal  angular  speed  (degree)'); 
meta  AVAZt 
/ipause 
prtsc 

%************************  Radar  Tracking  Parameter  output  ************ 
!del  tret. met 
plot(t.trel) .grid; 

title( 'RADAR  TRACKING  ELEVATION  ANGLE  degree'); 
xlabel( 1  In-flight  time  (sec)'); 
ylabel( 'Elevation  Angle  (degree)'); 

•meta  thel 
%  pause 
%  prtsc 

! del  trazl. met 
plot( t ,traz) .grid; 

title( 'RADAR  TRACKING  AZIMUTHAL  PARAMETER'); 
xlabel( ' In-flight  time  (sec)'); 
ylabel( ’ Azimuthal  angle  (degree)'); 
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meta  trazl 
»pause 
%  prtsc 

! del  AVAE.  met 
plot(avell , ) .grid; 

title( 'RADAR  TRACKING  ELEVATION  ANGULAR  VELOCITY  '); 
x 1 ab e 1 ( | MEASUREMENT  BASE  '); 

ylabel( 'Elevation  Angular  velocity  (degree/sec) 1 ); 
meta  AVAE 
X  pause 
X  prtsc 

! del  AVAZ. met 
plot(avazl , ) .grid; 

title( 'RADAR  TRACKING  AZIMUTHAL  ANGULAR  VELOCITY'); 
xlabel( 'MEASUREMENT  BASE'); 
ylabel( 'AZIMUTH  ANGULAR  (degree)’); 
meta  AVAZ 
*pause 
%prtsc 

EXETIMEl=etime( clock ,t2) 

%*■****■***■:>**************  graph  ************************************* 
!del  TVD.  met 
plot(t.fdl) .grid; 
tit le( 'DOPPLER  FREQUENCY'); 
xlabel( ' Flight  time  (sec)'); 
ylabel( ' Doppler  Frequency  (Hz)') 
meta  TVD 
X  pause 
prtsc 

! del  treVD. met 

plot( trel .Vdopl ) .grid; 

title( 'DOPPLER  Velocity'); 

xlabel( 'Radar  tracking  elevation  angle  (degree)'); 
ylabel( 'Doppler  Velocity(m) ' ) 
meta  trehVD 
%pause 
X  prtsc 

!del  trazVD.  met 

plot( traz , Vdopl ) , grid; 

title( 'DOPPLER  Velocity'); 

xlabel( ’ Radar  Tracking  azimuthal  Angle  (degree)'); 
ylabel( ' Doppler  Velocity(m) ' ) 
meta  trazVD 
%pause 
%prtsc 

!  del  RVDl.met 
plot(SRR, Vdopl) .grid; 
tit le( 'DOPPLER  Velocity'); 
xlabel( ' Slant  Range(m)'); 
ylabel( ' Doppler  Velocity(m) ' ) 
meta  RVD1 
%pause 
prtsc 

! del  tVDl.met 

plot(t .Vdopl) .grid; 

tit le( 'DOPPLER  Velocity'); 
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xlabel( ' In-flight  time(sec)') 
ylabel( 'Doppler  Velocity(m) ' ) 
meta  tVDl 
prtsc 

%************************  Radar  Tracking  Parameter  output  ************ 
! del  tret. met 
plot(t .trel) .grid; 

tit le( 'RADAR  TRACKING  ELEVATION  ANGLE  degree'); 
xlabel( ' In-flight  time  (sec)'); 
ylabel( ' Elevation  Angle  (degree)'); 
meta  thel 
%  pause 
prtsc 

! del  trazl. met 
plot(t .traz) .grid; 

title( fRADAR  TRACKING  AZIMUTHAL  PARAMETER'); 
xlabel( ' In-flight  time  (sec)'); 
ylabel( ’ Azimuthal  angle  (degree)'); 
meta  trazl 
«pause 
prtsc 

Vr  Vr  -.Y  -.V  -,V  tV  Vr  Vr  Vr  Ve  -jY  -iV  rr  Vr  i‘r  ->V  - V  ->Y  -jY  Vr  Vc  ■jV  Vr  Vr  Vr  Vr  Vr  *  Vc  Vr  Vr  Vc  VoY  -jY  tY  Vr  Vr  tWt  Vr  Vc  -jV  iY  Vr  *  Vr  ^Y  -ft  Vr  *  iY  *  Vr  Vr  *  *  iY  -jY 


%*******  RCS  and  Power  requirement  Estimation*****' 
load  50010 


* 


t3=clock 

Bvl=l. 2*max(fdl); 
Bv2=l. 2*22031; 
Bv3=l. 2*12721; 


%  video  Bandwidth  for  1st  filter 
%  video  Bandwidth  for  2nd  filter 
%  video  Bandwidth  for  3rd  filter 


%****************  radar  parameters  calculation  ************************ 


Loss=10**0. 24; 
KT=4. 021E-21; 
SN=10; 

Lamda=3E8/10525E6; 
Fn= 10**0.  2; 

G= 10**4.  43; 


%  total  loss  2. 4dB 
%  KT=K*T 

%  minimum  signal  to  noise  ratio  lOdB 
%  wavelength  meter  10525MHz 
%  noise  figure  2dB 
%  TX  or  RX  antenna  gain  44. 3dB 


const=(  ( 4*pi)**3*KT*Fn*Loss*SN)/(  ( G*Lamda )2 ) ; 

O/  >'■  j#  •* #  kV «*»  .v  »'«  j«  j at  hL 

/© 

%  Radar  cross  section  A  power  requirement  estimations 
for  n=l: N 

if  max(Height)=Height(n)  %  Summit  point  data 

ttop=t(n);  %  time  at  summit 

Ltop=SRR(n);  %  range  at  summit 

end 
end 

%*******************  switch  1  power  requirement  **************** 
Bcon=Bvl*const; 
for  n=l: 1092 
SRRa(n)=SRR(n); 
ta( n)=t( n) ; 
uDl=DR( : ,n); 
h=HeadingR( :  ,n); 
ct=h'*uDl; 

rcs=rcsa( ct);  %  res  function 

ctl(n)=ct; 
rcsl(n)=rcs; 
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thetal(n)=acos(ct)*180/pi; 

Ll=SRR(n); 

PtO=Ll**4*Bcon/rcs; 

Ptl(n)=PtO;  %  power  requied  in  free  space 

if  Ll=Ltop  %  required  power  at  summit 

ptop=PtO; 
end 

******************  *  switch  2  power  requirement  *  *******  *  *  ****** 
Bcon=Bv2*const; 
for  n=1093:  1630 
SRRa(n)=SRR(n); 
ta(n)=t(n); 
uDl=DR(: ,n); 
h=HeadingR( : ,n); 
ct=h'*uDl; 


rcs=rcsa( ct);  %  res  function 

ct l(n)=ct; 
rcsl(n)=rcs; 

thetal(n)=acos(ct)* 180/pi; 

Ll=SRR(n); 

PtO=Ll**4*Bcon/rcs; 

Ptl(n)=PtO; 
if  Ll=Ltop 
ptop=PtO; 
end 

switchl=max(Ptl); 
end  %  end  for  loop 

Ick'frit&'irsTiT’k'fri't'sriT'k'frk'i'rlrisltltlck  3rd  filter  ****%>*yrVf'>V,**^r**T,r*,/rVr**Vr**Vf 


%  power  requied  in  free  space 
%  required  power  at  summit 


Bcon=Bv3*const; 
for  n=1631:  1898 
SRRa(n)=SRR(n); 
ta(n)=t(n); 
uDl=DR( :  ,n); 
h=HeadingR( :  ,n); 
ct=h'*uDl; 


rcs=rcsa( ct);  % 

ctl(n)=ct; 
rcsl(n)=rcs; 

thetal(n)-acos( ct)* 180/pi; 
L1=SRR( n) ; 

Pt0=Ll**4*Bcon/rcs; 

Ptl(n)=PtO;  % 

if  Ll=Ltop  % 

ptop=PtO; 
end 

switch2=max( Pt 1 ) ; 
end  %  end  for  loop 
*£*********************  output 
! del  wwl 
diary  wwl 


res  function 


power  requied  in  free  space 
required  power  at  summit 


'k'k'jtfefrlt'frfrleirklrit'trk'klc'/t'ftirtrk'i'rfrlrkirk'kick'fck 


N=N 


f light time=max( t) 
stoptime=max( ta) 
stoprange=max( SRRa) 
stopfd=fdl( 1898) 
timetop=ttop 
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rangetop=Ltop 
power top=ptop 
EXETIM£2=et ime( clock , t 3 ) 


TOTALTIME=EXETIME 1  +  EXETIME2 


diary 

%****>■ r******************  graph  radar  cross  section  ***************** 
!del  rcsll.met 
plot(SRRa,rcsl , '  +  '  ) .grid; 

ylabel( ' Radar  cross  section  square  meters'); 
xlabel( ' Slant  range'); 
meta  rcsll 
«prtsc 
%  pause 
! del  rcsth. met 
plot( theta 1 , res  1 , '  +  ' )  , grid; 
ylabel(  '  Radar  cross  section  sejuare  meters'); 
xlabel( ' Aspect  angle  (degree)  j; 
meta  rcsth 
%  pause 
%prtsc 

! del  rasp,  met 
plot(ta,thetal , '+' ) .grid; 

ylabel( 'Radar  cross  section  sejuare  meters'); 

xlabel( ' Aspect  angle  (degree)  ); 

meta  rasp 

%pause 

%prtsc 

!del  rasl.met 

plot(ctl , '  +  ' ) .grid; 

xlabel( ' in-flight  time  (sec)'); 

ylabel( ' Aspect  angle  (degree)'); 

meta  rasl 

%prtsc 

°opause 

VrVr  *V?VoWr •frirfririrft'fc  Vr  5  fly  0  Variable  ^^*^ViV->yf***VrVr'>V-;V-A',5Vyf'jV,3V,j,rVrVr'>V'>vyrTV-5VVr'>V7V'3V?V-5V-V'sV 


%! del  Pw50010.  mat 

%save  Pw50010  ptlmax  ctl  t  fdl  rcsl  LI 
%**********************  printout  power 


N  thetal  Ptl 
requirement  *• 


SRR  N 


! del  powerl. met 
plot( SRRa , Pt 1 , ' +' ) .grid; 
title( 'RADAR  POWER  REQUIREMENT  '); 
xlabel( ' Slant  range  (m)  '); 
ylabel( ' Power  requirement  (Watts)’); 
meta  powerl 
prtsc 
%pause 

**********  RADAR  CROSS  SECTION  CALCULATION  ******************* 
clear 

LOAD  50010 
for  n=l: N 

uDl=DR(: ,n);  %  radar  to  target  unit  vector 

h=HeadingR( : ,n);  %unit  heading  vector  of  projectile 

ct=h'*uDl;  %  antenna  beam  angle 


^function  [ rcs.ct) =RCS(ct) 

%  constant  declaration 

b=0. 155/2;  %  shell  base  radius 
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1=0.5236;  %  shell  cone(nuzzle)  height 

i=0.4144;  %  shell  cylinder  height 

ta=b/l;  %  =atan(b/l)  i. e.  half  shell  cone  angle 

ta2=ta*ta; 

ca2=l/( l+ta2); 

ca=sqrt( ca2); 

cx=abs(ct); 

cx2=cx*cx; 

sx2=( l-cx2); 

sx=sqrt( sx2 ) ; 

rcsc=2*i*b,vsx; 

b2=b*b; 

rcsb=pi*b2*cx; 
if  cx>=ca 
rcs=rcsc+rcsb; 
else 

rt=ta*cx/sx; 

rt2=rt*rt; 

rcs2=b2';-(  asin(  rt)+rt*sqrt(  l-rt2) )*cx; 
res 1=( ( b*i) / ( sx2) )*( 1 -cx2/ca2 )**( 3/2 ) ; 
rcs=rcs l+rcs2+rcsc+rcsb/2; 

end  %  radar  cross  section  case  2 

rcs0(n)=rcs; 
ctl(n)=ct; 

thetal(n)=acos( ct)* 180/pi; 
end 

!del  rcsll 

%  plot( SKR, rcsO , 1 +  ' ) .grid; 

%ylabel( ' Radar  cross  section  square  meters'); 

»xlabel( ' Slant  range'); 

%meta  rcsll 

“oprtsc 

”0pause 
! del  rcsth 

plotf the tal , res 0, '+'), grid; 

ylabel( ' Radar  cross  section  square  meters'); 
xlabel( ' Aspect  angle  (degree)  ); 
meta  rcsth 
?*,pause 
! del  rasp 

plot(t , the tal , '+' ) ,grid; 

ylabel( 'Radar  cross  section  square  meters'); 

xlabel( ' Aspect  angle  (degree) r); 

meta  rasp 

%pause 

!del  rasl 

plot(ctl , '+’ ) .grid; 
xlabel( ' in-flight  time  (sec)'); 
ylabel( ' Aspect  angle  (degree)'); 
meta  rasl 

^******!’ov****vr*  end  of  program  ********************************* 
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APPENDIX  B.  CROSS-SECTION  ANALYSIS 

The  radar  cross-section  has  been  discussed  in  Chapter  4.  The  derivation  of  the  RCS 
equations  is  given  below.  The  radar  main  beam  direction  is  pointing  in  the  direction 
D.  The  projectile  axis  (or  heading)  is  pointing  along  h.  Accordingly,  the  coordinate 
system  (x,y.z)  can  be  chosen  so  that  x  =  D  and  the  y-z  plane  is  the  plane  on  which  the 
physical  cross-sectional  area  of  the  projectile  is  to  be  projected.  Another  coordinate 
system  (.v',y',z')  can  be  chosen  so  that  x'  =  h  and  the  projectile  is  a  body  of  revolution 
around  the  x'  axis.  These  two  coordinate  systems  can  be  set  up  jointly  as  follows  : 


'Q 

II 

(5.1) 

D  x  h 

\Dxh\ 

(5.2) 

'Q 

X 

"■u 

II 

<i* 

» 

(5.3) 

x'  —  h 

(5.4) 

K  1  1 

II 

(5.5) 

y'  =  ?'  x  v' 

(5.6) 

=  f  x  h 

(5.7) 

[(/i  x  D)h  -/)] 

(5.S) 

\Dxh 


I  '>  ~  rR  I 

Note  that  the  coordinate  (x.y.z)  are  used  in  this  Appendix  only  and  are  not  related 
to  other  parts  of  the  thesis. 

The  suface  of  the  projectile  is  given  by  : 
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s7+z'2  =  A*i 


(B.  10) 


The  RCS  is  the  area  enclosed  by  projecting  of  this  surface  on  the  y-z  plane. 

A.  CASE  1 

For  |  D  •  h\  >  cos  a.  Figure  25  on  page  38,  the  conoid  is  hidded  in  the  shadow  re¬ 
gion.  The  projected  areas  can  be  easily  seen  to  be  : 
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Ax  =  3.14159 

(B.  11) 

3 

II 

• 

(5.12) 

A2  =  lab 

(5.13) 

a2  =  A:  !  D  •  h  j 

(5.14) 

II 

■*T 

(5.15) 

o3  =  A3  j  D*  h  | 

(5.16) 

where  a  is  the  diameter  of  the  projectile  and  b  is  the  cylinder  height. 

area  in  : 

The  total  projection 

°i  =  °\  +  °2  +  ci 

(5.17) 

B.  CASE  2 

When  \D»h\<  cos  a.  the  geometry  of  the  target  is  shown  in  Figure  26  on  page 
39.  The  RCS  value  consists  of  four  different  parts:  A,  .  A:.  A,.  Aj .  Each  of  them  is  given 

below  : 

For  area  A,: 

y2  =  cV2 

(5.  IS) 

=  (■*'  -  x0')2  tan20asp 

(5.19) 

where 

c  =  tan  y. 

(5.20) 
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(c2  -  tan2easp).x'2  +  2x0’  ian26aspx'  -  x0'2  tan20flip  =  0 

(5.2 1) 

,  *o  tan20JV 

( tan2tfCJ/,  -  c2) 

(5.22) 

xc'c2  tan  easp 
( tan 26asp  -  c2) 

(5.23) 

where  xc',yc'  is  the  maximum  point  of  the  ellipse. 

c2 

x'  ~  tan  0asp 

(5.24) 

A*'*/)  =  S'  tan  6asf.  +  x' 

(5.25) 

X  =  V/ 

(5.26) 

(c2.v'  +  v'  tan 

( lan:0asp  cV  2 

(5.2~) 

At  a  point  which  is  p,  .  x’  =  L.  then: 

y'  tan  6asp  +  Lc 2  =  0 

(5.2  S) 

or 

lx2 

tan  0Jsp 

(5.29) 

So  the  coordinate  at  pi  is  {-L  ,  ,  0  )  on  the  z‘  axis, 

tan  0„, 

a,  * 

+2 -5  -  c2l: 

tan  8gSf 

(5.30) 

or 

cU  tan20aip  —  r2),  ;2 
tan  0asf, 

(5.31) 
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Hence 


^  t  ,  1  +  C  J  ;2 

A\  =  *L{ - T~ - ) 

tan  6 aSp 


=  L  c(  tan  9 asp  +  c  ) 


( tan"#™  -  o' ;2 


tan'tf 


asp 


Let  .V  be  the  normal  unit  vector  of  surface 


,V  = 


D  cos* a  —  /i  cos  9„ 


N  sin:#flV  cos4  a  ■+■  sin's  co$~9asp 


then  the  projection  area  for  .-1,  is: 


L 1  tan  z(  tan '9a.r  +  tan4*)1  2(  tan ~9a,r  -  tan's)1  *(#>  •  .V) 

tan'#™ 

u  » 


^  1 


1  -  2/ 


Oy  = 


nJ  COS*fla:r  ,  , 

aL  (1 - -  ~~~  )  " 


sin‘#av  cos's 


(5.32) 

(5.33) 


(5.34) 


(#.35) 


(5. 36l 


/l 2  is  part  of  a  semi-circle  (Figure  2*  on  page  39).  A  straight  forward  intergrauon  gives: 


A  2 


*>  ■  —  ] 
«'{  sin  ( 


tan  7 
;  tan  #c,. 


tan  s 
I  tan  9asr  | 


tan's 

tan'#™ 

u  t 


(5. 3") 


•  -i, 

a-  =  u  {  sin  ( 


tan  7 
I  tan  #av  I 


tan  s 
\  tan  9arr  1 


tan's 
tan  \Sf. 


}  I  0  •  A  I 


(5.3S) 


is  a  rectangle: 
A 3  =  2af> 
so 

o3  =  1  D  x  #  | 


(.#.39) 


(#.40) 
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where  a  is  the  radius  of  the  projectile  base,  and  b  is  the  height  of  the  cylindrical  part  of 
projectile. 

At  is  a  semicircle,  hence  : 


Aa  =  3.14159  — 


=  Aa\  D»h 


(B.A2) 


The  total  cross-section  a,  is  the  sum  of  upper  4  projection  areas. 


a,  =  a]  4-  n2  +  o3  +  <7  4 
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APPENDIX  C.  PROPAGATION  FACTORS 


Projectile 
Ground  Range 

Location 

Height 

10m 

Antenna 

15m 

Height 

20m 

30m 

50m 

impact  point  ,  , 

(m) 

O) 

(dB) 

(dB) 

(dB) 

(dB) 

(dtsj 

-482. 0000 

932. 4000 

0 

0 

0.  6000 

-0. 3000 

-0. 3000 

-382. 0000 

741.  2000 

o. 6000 

0. 6000 

0.  6000 

-0. 3000 

-0. 3000 

-333. 0000 

647. 9000 

-0. 7000 

-0. 4000 

0.  4000 

-0.  3000 

-0. 3000 

-234. 0000 

455. 2000 

-0. 3000 

-1.  1000 

1.  2000 

1. 2000 

-0. 2000 

-182. 0000 

355. 7000 

-0. 9000 

1.  0000 

1. 7000 

0.  7000 

1. 7000 

-152. 0000 

294. 8000 

1.  9000 

-3. 7000 

1.  2000 

-2. 6000 

1.  2000 

-148. 0000 

289. 2000 

2. 3000 

0.  5000 

-3. 5000 

1.  5000 

-3. 4000 

-145. 0000 

283. 7000 

1.  5000 

1. 0000 

0.  3000 

0 

0. 4000 

-142. 0000 

278. 1000 

-0. 7000 

1.  3000 

-1. 9000 

0. 3000 

-0. 7000 

-140. 0000 

272. 6000 

-3. 8000 

-2. 4000 

-1.  0000 

-1. 9000 

-2. 1000 

-137. 0000 

267. 0000 

-3. 0000 

-3. 0000 

-3. 0000 

-3. 2000 

-0. 9000 

-132. 0000 

261. 5000 

0 

1. 4000 

2. 3000 

2. 6000 

2. 3000 

-131. 0000 

256. 0000 

2. 3000 

2. 9000 

1. 5000 

-2. 4000 

-0. 3000 

-129. 0000 

250. 4000 

3. 2000 

0. 7000 

-4. 1000 

1.  3000 

-2. 1000 

-125. 0000 

244. 8000 

2. 2000 

-5. 1000 

0. 2000 

0. 2000 

2. 2000 

-123. 0000 

239. 3000 

-0. 8000 

-1.  6000 

3. 1000 

-1.  4000 

-0. 7000 

-120. 0000 

233. 7000 

-5. 0000 

2. 9000 

-4. 0000 

-5. 0000 

1. 6000 

-117. 0000 

228. 2000 

-4. 2000 

2. 9000 

-4. 0000 

-5. 0000 

1.  6000 

-114. 0000 

222. 6000 

-0. 1000 

-1.  0000 

2.  9000 

3. 2000 

0.  7000 

-111. 0000 

217. 1000 

3. 0000 

-6. 5000 

2. 7000 

-0. 3000 

-0. 3000 

-108. 0000 

211.  5000 

3. 7000 

0.  7000 

-5. 3000 

0. 8000 

-0. 3000 

-105. 0000 

206. 0000 

2. 7000 

3. 7000 

-0. 9000 

1. 8000 

2.  7000 

-103. 0000 

200. 4000 

-0. 3000 

2. 7000 

3. 8000 

-2. 4000 

-2. 4000 

-100. 0000 

194. 9000 

-6. 0000 

-2. 2000 

0. 6000 

3. 0000 

-0. 5000 

-97. 0000 

169. 3000 

-5.  3000 

-5. 3000 

-5. 3000 

-5. 0000 

3. 0000 

-94. 0000 

183. 8000 

-14. 7000 

-4.  7000  - 

•12. 5000 

3. 0000 

-11. 0000 

-91. 0000 

178. 2000 

3. 4000 

4. 4000 

3. 2000 

1.  4000 

-1.  7000 

-86. 0000 

172. 6000 

2. 2000 

4.  2000 

-3.  5000 

-0. 9000 

3.  1000 

-85. 0000 

167. 1000 

3. 5000 

-5. 5000 

-1.  5000 

3.  5000 

-0.  5000 

-83.  0000 

161. 5000 

0. 5000 

-3. 5000 

4. 2000 

-7. 5000 

-4. 8000 

-80. 0000 

156. 0000 

-6. 4000 

3.  4000 

2.  3000 

4. 1000 

2. 8000 

-77. 0000 

150. 4000 

4. 3000 

7. 1000  • 

-10. 0000 

-3. 0000 

3. 0000 

-74. 0000 

144. 8000 

0.  7000 

0. 9000 

2. 3000 

1.  9000 

-6. 9000 

-71.  0000 

139. 3000 

3.  7000 

-9. 9000 

4. 4000 

3. 0000 

3. 0000 

-68. 0000 

i33.  700u 

4.  6000 

-0. 4000 

-1.  4000 

-4. 4000 

3. 2000 

-66. 0000 

128.  1000 

3.  8000 

4. 4000 

-3. 3000 

4. 4000 

-6. 4000 

-63. 0000 

122.  5000 

1.  1000 

4.  1000 

4. 1000 

-4. 7000 

1.  9000 

-60. 0000 

117. 0000 

-6.  7000 

0.  1000 

3. 6000 

2. 5000 

4.  0000 

-57.  0000 

114.  4000 

-13.  4 000 

-4. 4000 

0.  2000 

4.  6000 

4.  6000 

-54. 0000 

105. 9000 

-0.  2000 

0. 2000 

0. 2000 

4.  1000 

2.  7000 

-52. 0000 

100.  3000 

4.  0000 

4. 8000 

4. 7000 

4.  7000 

3.  1000 

-49. 0000 

94.  7000 

5.  0000 

4, 0000 

2. 4000 

-3. 9000 

-6. 9000 

-46. 0000 

89. 2000 

4. 5000 

0.  2000 

-14. 4000 

2. 5000 

3. 6000 

-43. 0000 

83.  6000 

2. 8000 

-14. 5000 

1.  5000 

3.  5000 

1.  5000 

-40. 0000 

78. 0000 

1. 4000 

0.  2000 

4.  ;ooo 

-10.  2000 

-2. 2000 
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*37.0000  72.5000  -12.4000  4. 

-34. 0000  66. 9000  -4. 8000  4. 

-31.0000  61.3000  1.4000  1. 

-29.0000  55.7000  3.7000  -7. 

-26.0000  50.2000  4.7000  -4. 

-23.0000  44.6000  3.5000  2. 

-20.0000  39.0000  2.7000  3. 

-17.0000  33.4000  1.8000  2. 

-14. 0000  27. 9000  1. 1000  2. 

-11.0000  22.3000  -0.1000  0. 

-9.0000  16.7000  1.0000  -0. 

-6.0000  11.1000  -0.6000  -0. 

-3.0000  5.5000  4.7000  3. 

0  0  -2. 6000  -4. 


2000 

16. 0000 

4. 6000 

4.  6000 

7000 

-12. 6000 

0. 9000 

-3.  8000 

4000 

1.  4000 

-3. 1000 

2.  9000 

0000 

4.  3000 

4. 5000 

2.  6000 

6000 

2.  3000 

-1. 9000 

-3.  3000 

3000 

-6.  4000 

-1.  4000 

4.  5000 

7000 

-2.  1000 

4. 4000 

-11.  1000 

4000 

2.  4000 

0 

4.  5000 

4000 

3.  7000 

-1.  7000 

-9.  0000 

7000 

2.  2000 

3.  6000 

-10.  7000 

3000 

0. 3000 

3. 3000 

4.  1000 

9000 

-0. 2000 

-0. 2000 

4.  1000 

6000 

2. 7000 

2. 4000 

2.  8000 

8000 

-5. 0000 

-4. 6000 

-3.  2000 

APPENDIX  D.  MISSED  DETECTION  AND  FALSE  ALARM 


clear 

MN=1000; 

Npeak=MN/2; 
dtmax=0. 08650; 

Bvl=63900; 

FNnoise=Bvl*dtmax/Npeak; 

SNR0=10.  0; 
step=0. 5; 

Nf=5; 

Mf=29; 

thshstep=0.  005; 
thsh0=0.  355; 
for  m=l:  Mf 

thsh=m‘vthshstep+thshO; 
thresh( m)=thsh; 

FN£noise(  m)=FNnoise*2*asin(  sqrt(  thsh)  )/pi; 
end 

x=SNR0-step; 
for  n=l: Nf 
x=x+step; 

SNR( n)=x; 
sgnr=10**(0.  l*x); 

%  transfer  the  signal-to-noise  parameter  to  be  used  in  fdetect.m 
snr=sgnr 

!  del  hm  tempsnr.  mat 
save  hm  tempsnr  snr 

for  m=l: Mf 
thsh=thresh(m); 
pfan=exp( -thsh-sgnr); 

Pfa(n,m)=pfan; 

FNerrorC  n ,m)=pf an*FNEnoise(m) ; 
end 

Pmissfn , 1 )=Pmissdct( sgnr ,0 , thresh ( 1) ); 
for  m=2: Mf 

Pmissd=Pmissdct( sgnr ,thresh(m-l) ,thresh(m) ); 

Pmiss(n,m)=Pmiss(n,m-l)+Fmissd; 

end 

end 

! del  hm  pfaulse.mat 

save  hm  pfaulse  Nf  Mf  thresh  FNEnoise  SNR  Pmiss  FNerror  Pfa 

! del  hm  pfaulse 
diary  hm  pfaulse 

disp( ’  Threshold  Missed  Noise  Faulse’)^ 

disp^’  Detection  Error  Alarm  ') 

for  n= 1 : N f 
SN_ratio=SNR(n) 

(thresh;  Pmiss(n,:);  FNerror( n , : ) ;  Pfa(n,:))'; 
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end 

diary 

pause 

!del  hm  pfaulse. met 

plot  (thresh,  Pmiss( 1 ' ,  thresh,  FNerror( 1,:  ) , 'o' ) .grid; 

title( ' Signal-to-Noise  Ratio  =  10.0  db') 

xlabel( ' Threshold-to-Signal  Power  Ratio'); 

ylabel( ' Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

plot  (thresh,  Pmiss(2,: ) , '+' ,  thresh,  FNerror( 2 , : ) , ' o' ) , grid; 

title( ' Signal-to-Noise  Ratio  =  10.5  db' ) 

xlabel( ' Threshold-to-Signal  Power  Ratio'); 

ylabel( ' Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

plot  (thresh,  Pmiss( 3 ,:),'  +  ' ,  thresh,  FNerror( 3 , :  ) , ' o' ) , grid; 

title( ' Signai-to-Noise  Ratio  =  11.0  db') 

xlabel( ’ Threshold-to-Signal  Power  Ratio'); 

ylabel( ' Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

plot  (thresh,  Pmiss(4,: ) , '+' ,  thresh,  FNerror( 4 , :  ) , ' o' ) .grid; 

title( ' Signal-to-Noise  Ratio  =  11.5  db') 

xlabel( ' Threshold-to-Signal  Power  Ratio'); 

ylabel( ' Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

plot  (thresh,  Pmiss( 5 ,:),'+' ,  thresh,  FNerror( 5 , : ) , ' o' ) .grid; 

title( ' Signal-to-Noise  Ratio  =  12.0  db') 

xlabel( ' Threshold-to-Signal  Power  Ratio'); 

ylabel( ' Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 


********* 


-1.  J.  J. 


****************** 


*************************************** 


function  pmsd  =  Pmissdct( snr ,  thsha,  thshb) 

%  The  parameter  snr  for  the  singal-to-noise  ratio  (NOT  in  db)  should 
%  have  been  saved  to  hm  tempsnr. mat. 

%!del  hm  tempsnr. mat 
%save  hm  tempsnr  snr 

pmsd=quad( ' fdetect' , thsha , thshb , 0. 0001); 


*********** 


** 


’  * * *  * * ****** * * * ' 


******  *  *  *  *  **  **  ******  *  * 


%  The  parameter  sgnr  has  to  be  specified  in  the  program, 
function  pdens  =  fdetect(v) 
load  hm  tempsnr 
z iv=2*snr*s  q  r t ( v ) * i ; 

LIoz=log( BESSELN(0,ziv) ); 

Lpdens=LIoz+log( snr ) -snr*( v+1 ) ; 


76 


pdens=exp( Lpdens ) ; 

end  of  program 

%*****************  output  data  *******************Vr****Vr*********** * 
Threshold  Missed  Noise  Faulse 

Detection  Error  Alarm 


S/N  =  10.  0  db 


0. 3600 

0. 0266 

0. 1237 

0. 0273 

0. 3650 

0. 0278 

0. 1187 

0.  0260 

0. 37C0 

0. 0290 

0. 1138 

0. 0247 

0. 3750 

0. 0303 

0.  1091 

0.  0235 

0. 3800 

0. 0316 

0. 1046 

0.  0224 

0.  3850 

0. 0329 

0. 1002 

0. 0213 

0.  3900 

0.  0343 

0.  0961 

0.  0202 

0.  3950 

0. 0357 

0. 0921 

0.  0193 

0. 4000 

0.  0371 

0.  0883 

0.  0183 

0.  4050 

0. 0386 

0. 0846 

0.  0174 

0.  4100 

0.  0401 

0.  0810 

0.  0166 

0. 4150 

0. 0417 

0. 0777 

0.  0158 

0. 4200 

0. 0433 

0.  0744 

0. 0150 

0.  4250 

0. 0449 

0. 0713 

0.  0143 

0. 4300 

0. 0466 

0.  0663 

0.  0136 

0. 4350 

0. 0483 

0.  0654 

0.  0129 

0.  4400 

0. 0501 

0.  0627 

0. 0123 

0. 4450 

0.  0519 

0.  0600 

0. 0117 

0. 4500 

0. 0537 

0. 0575 

0.  0111 

0. 4550 

0. 0556 

0.  0551 

0.  0106 

0.  4600 

0. 0575 

0. 0527 

0. 0101 

0. 4650 

0. 0595 

0. 0505 

0. 0096 

0.  4700 

0.  0615 

0. 0484 

0. 0091 

0. 4750 

0. 0636 

0. 0463 

0. 0087 

0.  4800 

0. 0656 

0. 0443 

0. 0082 

0. 4850 

0.  0678 

0. 0424 

0.  0078 

0. 4900 

0. 0699 

0.  0406 

0. 0074 

0. 4950 

0. 0721 

0. 0389 

0.  0071 

0.  5000 

0. 0744 

0.  0372 

0. 0067 

S/N  =  10. 5 

db 

0.  3600 

0. 0211 

0. 0798 

0. 0176 

0. 3650 

0.  0221 

0.  0760 

0. 0166 

0. 3700 

0. 0232 

0. 0724 

0.  0157. 

0.  3750 

0. 0243 

0. 0690 

0. 0149 

0.  3800 

0. 0255 

0. 0658 

0. 0141 

0.  3850 

0. 0267 

0.  0627 

0.  0133 

0. 3900 

0. 0279 

0.  0597 

0.  0126 

0. 3950 

0.  0291 

0. 0569 

0. 0119 

0.  4000 

0. 0304 

0. 0542 

0. 0112 

0. 4050 

0. 0318 

0. 0516 

C. 0106 

0. 4100 

0.  0332 

0.  0491 

0.  0100 

0. 4150 

0. 0346 

0. 0468 

0. 0095 

0.  4200 

0.  0360 

0. 0446 

0.  0090 

0. 4250 

0. 0375 

0. 0424 

0. 0085 

0.  4300 

0.  0391 

0. 0404 

0. 0080 

0.  4350 

0. 0406 

0. 0385 

0. 0076 
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0.4400  0.0423 

0. 0366 

0. 0072 

0.4450  0.0439 

0. 0349 

0. 0068 

0.4500  0.0456 

0. 0332 

0. 0064 

0.4550  0.0474 

0.  0316 

0. 0061 

0.4600  0.0492 

0.  0301 

0. 0057 

0.4650  0.0510 

0. 0286 

0. 0054 

0.4700  0.0529 

0. 0272 

0. 0051 

0.4750  0.0548 

0. 0259 

0. 0048 

0.4800  0.0567 

0.  0247 

0. 0046 

0.4850  0.0588 

0. 0235 

0. 0043 

0.4900  0.0608 

0.  0223 

0. 0041 

0.4950  0.0629 

0. 0213 

0. 0039 

0.  5C00  0.  0650 

0. 0202 

0. 0037 

S/N  =  11.  0  db 

0.3600  0.0163 

0. 0487 

0.  0108 

0.3650  0.0172 

0.  0461 

0. 0101 

0.3700  0.0181 

0. 0436 

0. 0095 

0.3750  0.0191 

0.  0413 

0. 0089 

0.3800  0.0201 

0.  0391 

0. 0084 

0.3850  0.0211 

0. 0370 

0. 0079 

0.  3900  0.  0222 

0. 0350 

0. 0074 

0. 3950  0. 0233 

0. 0331 

0. 0069 

0.4000  0.0244 

0. 0313 

0. 0065 

0.4050  0.0256 

0. 0296 

0. 0061 

0.4100  0.0268 

0. 0280 

0. 0057 

0.4150  0.0281 

0. 0265 

0. 0054 

0.4200  0.0294 

0.  0251 

0. 0051 

0.4250  0.0307 

0. 0237 

0. 0047 

0.4300  0.0321 

0. 0224 

0.  0045 

0.4350  0.0335 

0.  0212 

0.  0042 

0.4400  0.0350 

0.  0201 

0. 0039 

0.4450  0.0365 

0.  0190 

0. 0037 

0.4500  0.0360 

0.  0179 

0.  0035 

0.4550  0.0396 

0. 0170 

0. 0033 

0.4600  0.0413 

0.  0160 

0.  0031 

0.4650  0.0430 

0.  0151 

0. 0029 

0.4700  0.0447 

0.  0143 

0.  0027 

0.4750  0.0465 

0.  0135 

0. 0025 

0.4800  0.0483 

0.  0128 

0.  0024 

0.4850  0.0502 

0. 0121 

0.  0022 

0.4900  0.0521 

0. 0114 

0.  0021 

0.4950  0.0540 

0.  0108 

0. 0020 

0.5000  0.0561 

0.  0102 

0.  0018 

S/N  =  11. 5  db 

0.3600  0.0123 

0.  0280 

0.  0062 

0.3650  0.0130 

0.  0263 

0.  0058 

0.3700  0.0138 

0.  0247 

0.  0054 

0.3750  0.0146 

0.  0232 

0.  0050 

0.3800  0.0154 

0.  0218 

0.  0047 

0.3850  0.0163 

0.  0205 

0.  0043 

0.3900  0.0172 

0.  0192 

0.  0041 

0.3950  0.0182 

0.  0181 

0.  0038 
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0. 4000 

0.  0191 

0. 0169 

0. 0035 

0.4050 

0.  0201 

0. 0159 

0. 0033 

0.  4100 

0. 0212 

0. 0149 

0.  0031 

0.4150 

0. 0223 

0. 0140 

0. 0028 

0. 4200 

0. 0234 

0. 0132 

0. 0027 

0. 4250 

0. 0246 

0. 0123 

0. 0025 

0.  4300 

0.  0258 

0. 0116 

0. 0023 

0. 4350 

0. 0271 

0. 0109 

0. 0021 

0.  4400 

0. 0284 

0.  0102 

0.  0020 

0. 4450 

0. 0297 

0. 0096 

0. 0019 

0.  4500 

0.  0311 

0. 0090 

0. 0017 

0. 4550 

0. 0325 

0. 0084 

0. 0016 

0. 4600 

0. 0340 

0. 0079 

0. 0015 

0. 4650 

0. 0355 

0. 0074 

0. 0014 

0.  4700 

0. 0371 

0. 0070 

0. 0013 

0. 4750 

0. 0387 

0. 0065 

0. 0012 

0.  4800 

0. 0404 

0. 0061 

0. 0011 

0. 4850 

0.  0421 

0. 0057 

0. 0011 

0.  4900 

0.  0439 

0. 0054 

0.  0010 

0.  4950 

0. 0457 

0. 0050 

0. 0009 

0. 5000 

0. 0476 

0. 0047 

0. 0009 

S/N  =  12.  0  db 

0. 3600 

0.  0090 

0.  0151 

0. 0033 

0. 3650 

0. 0096 

0. 0140 

0.  0031 

0. 3700 

0.  0102 

0.  0131 

0. 0028 

0.  3750 

0.  0108 

0.  0122 

0. 0026 

0. 3800 

0.  0115 

0.  0113 

0. 0024 

0. 3850 

0. 0122 

0. 0105 

0. 0022 

0. 3900 

0.  0130 

0. 0098 

0.  0021 

0. 3950 

0. 0138 

0.  0091 

0. 0019 

0.  4000 

0.  0146 

0. 0085 

0. 0018 

0.  4050 

0.  0154 

0. 0079 

0. 0016 

0. 4100 

0. 0163 

0. 0074 

0. 0015 

0.  4150 

0. 0173 

0. 0069 

0. 0014 

0. 4200 

0.  0182 

0. 0064 

0.  0013 

0.  4250 

0.  0192 

0.  0059 

0.  0012 

0. 4300 

0. 0203 

0. 0055 

0. 0011 

0. 4350 

0. 0214 

0.  0051 

0.  0010 

0. 4400 

0. 0225 

0. 0048 

0. 0009 

0. 4450 

0.  0237 

0. 0044 

0.  0009 

0. 4500 

0. 0249 

0. 0041 

0. 0008 

0.  4550 

0. 0261 

0.  0038 

0. 0007 

0. 4600 

0. 0275 

0. 0036 

0.  0007 

0. 4650 

0.  0288 

0. 0033 

0.  0006 

0. 4700 

0. 0302 

0. 0031 

0.  0006 

0. 4750 

0. 0317 

0. 0029 

0. 0005 

0. 4800 

0. 0331 

0. 0027 

0. 0005 

0.4850 

0.  0347 

0. 0025 

0. 0005 

0. 4900 

0.  0363 

0. 0023 

0. 0004 

0. 4950 

0.  0379 

0. 0022 

0. 0004 

0. 5000 

0. 0396 

0.  0020 

0. 0004 

yjrklrick&Jeirirlcklrirlrlrlcirk'ir 

end  of  data 
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